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Abstract 


A  rotating,  cylindrical  combustion  chamber 
was  used  to  observe  the  effect  of  chamber  rotative  speed 
and  location  of  the  point  of  ignition  on  combustion  duration 
and  flame  trajectory  in  a  cylindrical  force  field.  Experi¬ 
mental  trajectories  of  the  flame  in  homogeneour  and  non- 
homogeneous  mixtures  were  compared  with  analytically 
obtained  ones.  The  study  lead  to  the  following  conclusions: 

1.  The  optimum  location  of  the  point  of  ignition  is  at  the 
periphery  of  the  combustion  chamber  for  both,  homogeneous 
and  nonhomogeneous ,  gas  mixtures  because  this  flame 
origin  produces  the  strongest  vortex. 

2.  Combustion  of  a  nonhomogeneous  gas  mixture  provides 
g r v  at e r  time  rates  of  energy  emission  than  that  of  a  homo¬ 
geneous  mixture.  Combustion  of  a  homogeneous  mixture, 
however,  produces  a  stronger  swirl  than  that  of  a  non¬ 
homogeneous  mixture  when  both  are  overrich. 

3.  At  the  beginning  of  the  combustion  the  experimentally 
observed  trajectories  of  the  flame  coincide  with  the  theo¬ 
retically  determined  trajectories,  but  in  the  later  portion 
of  the  combustion  they  do  not  because  the  drag  on  the 
inflamed  portion  of  the  mixture  exceeds  that  predicted  by 
theory. 
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I.  INTRODUCTION 

1.  1  General  Statement  of  the  Problem 

In  recent  years,  several  investigators  have  pointed  out 
that  in  a  vortex  type  combustion  chamber  the  inflamed  gas 
particles,  being  of  lower  density  than  the  surroundings,  have 
the  tendency  to  move  spirally  toward  the  center  of  the  chamber 
under  the  influence  of  the  centrifugal  force  field  of  the  air 
swirl.  If  we  are  to  enter  into  the  field  of  research  pertaining 
to  the  thermal  mixing  effect  of  gases,  then  it  is  necessary  to 
know  the  flow  pattern,  at  least  up  to  the  moment  at  which  com¬ 
bustion  starts,  since  all  further  particle  motion  will  be  governed 
by  the  existing  flow  patterns.  It  is  difficult  to  predict  the  flow 
pattern  in  an  actual  engine  because  the  moving  piston  produces 
transient  flows  that  result  in  random  turbulence.  Also,  there 
are  wall  effects  that  at  present  are  not  fully  understood. 

It  is  possible  to  eliminate  the  effect  of  the  combustion 
chamber  wall  by  rotating  the  chamber  at  constant  speed  and 
thus  produce  a  flow  pattern  which  approaches  that  of  perfect 
solid  body  rotation.  For  this  reason,  this  study  was  under¬ 
taken  using  a  rotating  cylindrical  combustion  chamber  in 
order  to  determine  the  trajectory  of  the  flame.  It  was  con¬ 
sidered  of  value  to  study  the  effect  of  chamber  speed  and 
spark  plug  position  on  the  combustion  time  in  both  homo¬ 
geneous  and  nonhomogeneoue  gas  mixtures.  In  order  to 
observe  and  study  these  phenomena,  flame  photographs  were 
taken  with  a  high  speed  camera. 
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The  ultimate  object  of  the  study  was  to  compare  the  theo¬ 
retical  and  experimental  flame  trajectories,  that  is,  the  locus  of 
the  center  of  the  flame.  The  assumption  was  made  that  only  the 
pressure  force  due  to  the  centrifugal  action  and  the  drag  force 
on  the  flame  significantly  influence  its  path.  Magnetic  forces 
and  other  miscellaneous  effects  were  neglected  since  the  chamber 
was  not  operated  at  elevated  temperatures. 

1 .  2  Previous  Related  Studies 

F.  and  A.  Pischinger  (1)  treated  this  phenomenon  theo¬ 
retically  on  the  assumption  that  the  air  rotation  consists  of  a  free 
vortex.  They  came  to  the  conclusion  that  the  trace  of  the  flame 
is  a  logarithmic  spiral,  the  angle  of  which  depends  on  the  ratio 
of  burnt  gas  density  to  that  of  the  unburnt  gas,  and  also  that  the 
angle  is  independent  of  the  swirl  speed.  F.  Naga  >  and  H. 

Kakimoto  (2)  investigated  a  spiral  curl  of  flame  in  a  swirl  chamber 
and  discussed  the  reason  why  better  combustion  results  when  the 
fuel  was  injected  in  the  direction  of  the  air  swirl  rather  than 
against  it.  J.  S.  Meurer  (3)  emphasized  in  a  recent  paper  that 
the  thermal  mixing  effect  was  one  of  the  most  important  factors 
that  let  the  M-system  approach  the  ideal  combustion  process. 
Ei&feld  (4)  determined  the  trajectory  of  droplets  within  the  Stoke's 
regime  in  a  swirling  field.  The  problem  of  describing  the  tra¬ 
jectories  of  particles  has  also  been  solved  by  Kriebel  (5)  by  means 


Numbers  in  parentheses  designate  References  in 
Bibliography. 
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of  the  Laplace  transformation  under  the  assumptions  that  there  is 
a  density  difference  between  the  surroundings  and  the  particles 
and  that  the  only  factors  that  affect  the  trajectories  are  the 
pressure  forces  due  to  the  centrifugal  action  and  the  drag  forces* 
Calder  (6)  investigated  the  motion  of  a  heavy  spherical  particle 
in  the  rotating  cylinder.  He  assumed  that  the  viscous  drag  on  the 
particle  was  given  by  the  Stokes  resistance  formula,  and  that  the 
fluid  rotated  as  a  solid  body* 

The  flow  pattern  in  a  model  chamber  with  a  tangential 
inlet  passage  was  observed  and  measured  by  A.  W.  Hussmann 
and  F.  Kahoun  (7)  who  determined  that  the  flow  was  one  of  solid 
body  rotation  and  was  not  stable. 

1 .  3  Origin  of  the  Study 

A  few  theoretical  studies  have  been  published  relative  to 
the  combustion  of  a  uniform  mixture  in  a  static  vessel.  S.  D. 
Raezer  (8)  calculated  the  relationship  between  burning  velocity 
and  space  velocity  of  a  spherical  wave  in  a  fixed  and  closed 
spherical  chamber.  A.  S.  Campell  (9)  studied  the  required 
combustion  time  for  constant  volume  combustion  in  the  one- 
dimensional  case.  F.  Nagao  et  al  (10)  discussed  the  flame 
development  in  a  swirl  with  solid  body  rotation.  D.  G.  Udelson  (11) 
experimentally  studied  the  effect  of  a  single  sphere  of  burning 
vapors  resulting  from  a  drop  of  liquid  fuel  in  an  undisturbed 
air  stream,  but  he  did  not  consider  the  force  effects  due  to  an 
unbalanced  system,  nor  did  he  investigate  flame  development 
in  a  wheel-type  swirling  field. 
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Thermal  mixing  is  an  important  factor  in  the  improvement 
of  fuel  stratification  and  combustion.  The  effect  of  thermal  mixing 
on  combustion  is  known  to  be  related  to  the  swirling  motion  of  the 
field  (12).  After  ignition,  this  swirl  will  be  intensified  by  the  flame 
motion.  Therefore,  one  important  factor  in  investigating  the 
thermal  mixing  is  to  find  the  trajectory  of  the  flame  in  a  swirling 
field  after  ignition  has  occurred. 

In  this  investigation  the  effect  of  radial  unbalance  was 
inclqded  in  the  theoretical  determination  of  the  flame  trajectory 
for  both  fixed  and  rotational  coordinates  under  the  assumptions 
of  solid  body  rotation  and  negligible  wall  effect.  The  flame  tra¬ 
jectory  was  also  investigated  by  experimental  methods. 
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IL  THEOR.ETJCAL  CONSIDERATIONS 

2.  1  Statement  of  the  Problem 

In  this  investigation  an  analysis  has  been  made  of  the  flame 
trajectory  in  a  rotating  combustion  chamber  model  during  the 
combustion  process.  The  problem  is  very  difficult  to  analyze 
precisely  because  of  the  many  unknown  functions  of  time.  For 
example,  the  relationships  between  pressure  and  time  and 
between  the  radius  of  the  flame  and  time  are  important  factors 
in  finding  the  flame  trajectory,  but  these  relationships  are  very 
difficult  to  obtain  analytically  and  must  be  determined  experi¬ 
mentally.  Von  Elbe  (13),  using  the  concept  of  burning  fraction 
which  was  derived  by  A.  S.  Campell  (9)  has  analytically  pre¬ 
dicted  the  radius  of  a  spherical  flame  under  the  assumption  that 
the  flame  expands  adiabatically .  Von  Elbe's  (13)  equation  cannot 
be  used  for  the  present  problem  because  the  relationship  between 
pressure  and  time  is  an  unknown  function.  In  this  study  the  wall 
effect  was  neglected,  that  is,  the  derivation  that  will  be  pre¬ 
sented  holds  only  until  the  flame  front  reaches  the  cylinder  wall. 

Since  the  drag  coefficient  is  an  unknown  function  of  the 
Reynolds  number,  the  relationship  between  drag  coefficient  and 
Reynolds  number  was  approximated  in  the  present  analysis  by 
a  piecewise  logarithmic -linear  relation.  In  the  latter  portion 
of  the  combustion  time,  as  shown  in  Fig.  2.  1,  indistinct 
boundaries  between  the  flame  and  its  surroundings  occur.  The 
flame  core  consists  of  decomposition  products  of  acetylene  and 
air,  and  the  amorphous  portion  may  be  considered  to  be  a  kind 
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of  flame  diffused  from  the  flame  core.  Af  ter  the  flame  grows  to 
an  order  of  magnitude  comparable  to  that  of  the  combustion 
chamber,  drag  force  is  no  longer  proportional  to  the  square  of 
relative  velocity.  Therefore,  we  cannot  apply  the  piecewise 
logarithmic  linear  relationship  to  this  situation;  thus  the  assumed 
relationship  between  drag  coefficient  and  Reynolds  number  does 
not  hold  for  the  later  portion  of  the  combustion  period. 


Fig.  2.  1.  Schematic  diagram  of  the  growing  flame 
2.  2  As  sumptions 

In  the  following  simplified  analysis,  the  resulting  flame 
trajectory  of  an  acetylene -air  mixture  injected  into  a  rotating 
cylindrical  chamber  filled  with  air  and  ignited  by  a  spark  plug 
has  been  investigated  under  the  following  assumptions: 

1.  The  flame  is  two-dimensional;  in  other  words, 
combustion  flame  spreads  only  in  the  plane 
perpendicular  to  the  axis  of  the  combustion 
chamber. 

2.  The  radius  of  the  flame  is  constant  over  a  small 


interval  of  time. 
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3.  The  angular  velocity  of  the  chamber  is  constant. 

4.  The  boundary  line  between  the  flame  and  its 
surrounding  is  distinct. 

5.  The  flame  moves  through  an  acetylene -air  mixture 
that  is  homogeneous  and  rotates  as  a  solid  body. 

6.  Magnetic  forces  are  negligible. 

2.  3  Derivation 

The  equation  of  equilibrium  is  as  follows: 

Inertia  force  =  drag  force  +  radial  pressure  force 


F  =  F~  +  F 
a  D  p 


Pressure  force:  When  the  flame  size  is  small,  the 

pressure  gradient  in  the  radial  direction 
is  Sol  X •  When  the  flame  grows 
to  a  size  comparable  to  the  chamber 
size,  the  pressure  force  in  a  radial 
direction  must  be  obtained  by  a  different 
analysis. 


Fig.  2.  2.  Radial  pressure  force  diagram  for  large 
flame  size. 
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As  noted  in  Fig.  2.  2,  since  the  circumferential  and  axial 
pressure  gradients  are  negligible,  one  can  obtain  the  radial 
pressure  force  applied  to  the  flame  from  the  Navier -Stokes 


equation  in  cylindrical  coordinates. 

The  radial  pressure  gradient  is 


?=  T]£rluj2+C 

Letting  P  be  the  pressure  above  the  pressure  at  the 
geometrical  center,  then 

c  =  0 

Now  the  pressure  force  is  given  by 

FP=  $>T>cos&<]/\ 

where 

J/|—  Rhd.0 

Therefore 

f>  =  2/*^  fa.  rWRh  cose)<& 

___  JU  faJth  K 

~ 

The  flame  volume  is  given  by 

V  =. 

Therefore  the  radial  pressure  force  can  be  expressed  as 


f'f9  Ja. 
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Thus,  even  in  the  case  of  a  large  size  flame,  the  same 
equation  for  the  radial  pressure  force  can  be  used  as  in  the  case 
of  a  small  size  flame. 

The  drag  force  is  given  by 

/ 


Since  the  drag  force  and  the  pressure  force  are  directed 
inward  toward  the  geometrical  center,  the  equation  of  equilibrium 
c?nbe  rewritten  as 

F  =  -Fn  -F 
a  Dp 

The  mas 8  of  the  flame  is 

Or)j_=  Z 

Therefore,  the  equilibrium  equation  can  be  shown  to  be 


cU J«.R.  h  Cp 

S7CF  —~~f - 


Htel  - 
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Z  u)“ 


Dividing  both  sides  of  results  in  the  following  expression 
for  the  acceleration  of  the  flame: 


di-*-  xKfj"'  ff  * 


Eq.  2. 1 


Now 


fq_  CZ  a  r\ 

T-A0  for  the  acetylene-air  mixture  flame;  therefore, 
in  this  case,  Equation  2.  1  is  as  follows: 
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Fig.  2.3  Vector  diagram  for  relative  velocity  of 
the  flame. 

As  obtained  from  Fig.  2.  3,  the  relative  velocity  of  the 
flame  can  be  expressed  as 

^Y&l  =  ^  Eq.  2.3 

But 

Wxrt=M.(UiJf)xvK 

=  -  tcoXj-w^ 

=  'I  g  cO  Eq.  2.4 

Substitution  of  Equation  2.4  into  Equation  2.  3  results  in: 

\4,= 
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Substitution  of  this  relation  into  Equation  2.  2  results  in 
the  following  expression  for  the  acceleration  of  the  flame 


djL.—  (g.jggj)  -3S  uJ'Z  Eq.2.5 

d  tA  XR  ' 


If  one  changes  to  a  rotational  coordinate  system,  the 
tangential  velocity,  which  is  given  by  0  —  Vo  in  the  fixed 
coordinate  system,  will  be  modified  as  follows 

ds  iiw.  o|K  \ 

KJ fldt  x  wt  v 


—  +•  UJ*"J 

•  — 

=  _oj  xtt) 


Then  for  the  two  dimensional  case 


The  acceleration  for  the  fixed  coordinate  system  is  given 


‘  *  _  v  —  — 

=  J2«t>  x«j  )4w  X*) 

In  the  present  problem 

cJ  =  0 


12 


Therefore,  for  the  two  dimensional  case 

»  • 

&<,=  ^  <-  «U  -  <y]  CO2- 

If  one  breaks  the  complex  numbers  into  their  real  and 
imaginary  parts,  the  results  are: 

for  fixed  coordinates 


for  rotational  coordinates 


■S  =  K  +  l 


or 


therefore 


xr  \  „  id 


z  —  {"\  t ivoy)e 


Eq.  2.  6 


iuSt 


Eq.  2.  7 


Substitution  of  Equations  2.  6  and  2.  7  into  Equation  2.  5 
will  give: 


tu)l 


*)  £iUjt 


-301  u)xeia}^ 
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Finally,  there  results  in  the  rotational  coordinate  system 


•  *  • 


Eq.  2.  8 


In  Equation  2.  8  the  first  term  of  the  left  side  is  an  apparent 
acceleration  term  of  the  flame  center.  The  second  term  of  the 
left  side  is  the  Coriolis  force  component  due  to  rotation. 

The  boundary  conditions  for  Equation  2.8  are 

^=-0  at  t  —  o 


Since  the  axis  through  the  point  of  ignition  was  taken  as 
the  real  axis,  the  point  of  ignition  is  a  real  value.  Therefore, 
for  a  cylinder  of  2.  5  inches  diameter  with  the  point  of  ignition 
midway  between  the  periphery  and  the  center 

fte  ^  =  o.  6  j  -vtt)  'vj"  =.  q 

Separating  "J  into  its  real  and  imaginary  parts,  one  has 
for  the  boundary  conditions  in  the  rotational  coordinate  system: 

§.  =  nh  f «  o 

(B) 

=0.^-2 5^  £e=0  CLi  h— ‘O 


Recalling  Equation  2.  5  and  2.  6  in  order  to  determine  the 
boundary  conditions  in  the  fixed  coordinate  system 
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When  one,  therefore,  has: 

il  =  Ho 

1*  =  l*% 

Separating  and  °*J  0  into  their  real  and  imaginary 

parts,  one  has: 

tfo  +  i  j.— <»(*.  + *5.) 

Therefore,  the  boundary  conditions  for  fixed  coordinates 

are: 

^0=  t?  at  -t  =  o  j 

Xo=0,  j-0.«5W  «tt=°  J 

To  go  further  with  Equations  2.  5  and  2.  8  we  need  to  know 
the  drag  coefficient  and  since  it  is  related  to  the  Reynolds 
number,  Re,  it  is  necessary  to  find  the  Reynolds  number  which 
is  given  by  the  following  expression. 


because  Re  is  a  scalar  value. 


But 

14 1  =  x-iz60 
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Where  l//«  I  is  the  complex  conjugate  of 
Therefore,  one  has  ^ 

Using  this  expression  for  the  relative  velocity,  Re  is 
given  by: 

1 

f\*6.  — —  (  0(  ~f(  ^  in  fixed  coordinates. 

Now  in  rotational  coordinates: 

0 

*  X 

=  54<? 

Therefore,  the  Reynolds  number  in  rotational  coordinates 
is  given  by: 


It  should  be  noted  that  y^j  and  D  are  scalar 

values;  thus  these  parameters  are  invariants  in  the  coordinate 
transformation.  Since  the  relationship  between  drag  coefficient 
Cp  and  the  Reynolds  number  Re  for  a  cylindrical  object  is  not 
defined  by  a  simple  function,  this  relationship  was  approximated 
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by  a  series  of  logarithmic  linear  functions  as  shown  in  Table  4.  1 
and  Fig.  4.  16. 

If  Equation  2.  5  is  separated  into  its  real  and  imaginary 
parts,  one  has: 


-3 


Eq.  2.9 


Similarly,  for  Equation  2.  8: 


i= 


^&..rST7T 


5  =^20lS  - 

s - -Jig  i  -JS^’ 


Eq.  2.  10 


Equation  2.  9  and  2.  10  are  the  fundamental  differential 
equations  for  this  investigation;  from  these  equations  the  tangential 
velocity,  the  position  vector,  and  thus  the  trajectory  of  the  flame 
can  be  calculated  under  the  assumptions  given  in  Section  2.  2.  It 
should  be  noted  that  in  Equations  2.9  and  2.  10,  the  flame  radius  R 
as  a  function  of  time  must  be  calculated  from  measured  values  of 
the  instantaneous  combustion  pressure  P  and  from  its  final  value 


P  .  Also  Cp  must  be  found  from  Re  after  finding 
measured  values  of  Pe  and  P  and  combustion  time. 


M  I? 


using 
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Boundary  conditions  (B)  for  rotational  coordinates  apply 
to  Equation  2.  10. 

Boundary  conditions  (C)  for  fixed  coordinates  apply  to 
Equation  2.9. 

To  determine  the  radius  of  the  flame  for  use  in  Equations 
2.  9  and  2.  10,  the  method  of  B.  Lewis  and  G.  Von  Elbe  (13)  used 
to  derive  an  expression  for  the  radius  of  a  spherical  flame  was 
applied  to  a  cylindrical  flame.  The  derivation  shown  in  Appendix  I 
results  in  the  following  expression  for  the  radius  of  the  cylindrical 
flame: 


r =iu 


Eq.  2.  11 


Where,  as  shown  in  Appendix  I 

Pc  -po  Eq.  2.  12 

Pq  is  the  initial  pressure  of  the  mixture  before  ignition 

occurs. 

Equation  2.  12,  derived  from  energy  considerations  and 
conditions  in  the  burnt  gas  region,  is  valid  only  in  the  early 
stages  of  burning,  in  which  the  pressure  rise  is  small.  S.  W. 
Raezer  (8)  has  derived  the  final  pressure  thermodynamically 
for  combustion  at  constant  volume,  using  average  specific  heats 
and  dissociation  equilibria.  In  the  present  study,  P^  and  P  were 
measured  by  a  pressure  transducer  and  an  oscilloscope.  Thus 


it  is  possible  to  calculate  the  flame  radius  R  from  E 
using  measured  values  of  P. 

Since  the  surrounding  mixture  will  be  compr 
flame  expands  adiabatically,  the  temperature  and  d< 
mixture  will  also  change  adiabatically.  Also,  since 
of  the  mixture  is  a  function  of  the  temper ature,  it  w 
from  the  air  tables  (neglecting  the  pressure  effect), 
ships  between  radius  of  the  flame  and  time  and  betw 
(recall  the  Reynolds  number)  and  time  were  calcula 
Equation  Z.  11  and  the  air  tables.  The  change  of  rai 

P  f> 

and  of  WH1  be  described  in  detail  in  Cl 

After  the  above  factors  were  found,  Equatioi 
were  solved  by  numerical  methods  (Kutt  a -Simpson  ■ 
with  the  aid  of  the  University's  digital  computer,  ar 
The  derivation  of  the  numerical  analysis  for  these  € 
shown  in  Appendix  I;  the  computer  program  is  show 
II. 

2.  4  Hypothetical  Considerations 

In  order  to  understand  more  clearly  the  fore 
we  will  consider  two  simplified,  hypothetical  cases 
These  cases  are: 

(1)  The  only  force  acting  on  the  systerr 
pressure  force. 

(2)  The  only  force  acting  on  the  systerr 
drag  force. 
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One  considers  these  cases  in  order  to  determine  the 
relative  magnitude  of  the  pressure  force  and  drag  force  acting  on 
the  flame. 

Case  1 

Consider  the  case  in  which  the  only  force  acting  on  the 
system  is  the  pressure  force.  In  this  case,  the  equation  of  motion 
for  a  fixed  coordinate  system  is  given  by 

1*—  -3  ^ 

Its  solution  is 

c,  4  C2 


The  following  results  can  be  obtained  by  applying  the 
boundary  conditions  for  the  fixed  coordinate  system,  given  in  (C), 
to  ^  and  its  first  derivative: 


c, + Cz  —  o .  62  jr 
c  c  —  # 


C, 

$  -  1 0(^-0 


d 


*VT 

0.  CoS 


3  -  &***"•>* 
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If  one  rewrites  the  above  equations,  one  has: 


or 


Eq.2.  13 


(o.6*sy-  (o-ttsfirT 


For  a  rotational  coordinate  system  we  have: 

S  —  ^=—‘({3.C^\lTujh,Zos Lvt  -+ S^V^wtS^t^ 


Eq.2.  14 


^  \fy  ^■Si'Ml/Jwt.Ceseot  5  Co3\f3U)t. 


IS 


If  one  writes  Equation  2.  13  in  general  form,  the  result 


^  r 


/i2  (V/B-r 


Eo  2.  15 


where: 

A  =  position  of  the  spark  plug 

B  =  &/Sk 


Case  2 

Here  the  pressure  force  is  neglected  and  only  the  drag 
force  is  assumed  to  be  acting  on  the  flame.  In  this  case,  the 
resulting  equations  of  motion  of  the  flame  in  rectangular 


coordinates  are 
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j£ _ 

1=  "XR 


Eq.  2.  16 


Equations  2.  16  were  solved  on  the  I.B.M.  7074  computer. 
When  the  position  of  the  spark  plug  is  fixed  and  the  density- 
ratio  of  mixture  to  flame  is  increased,  Equation  2.  15  forms  a 
family  of  elliptical  curves.  When  — ^ - r  fro  ,  the  flame 

5a- 

motion  is  that  of  simple  harmonic  motion.  A  more  detailed 
discussion  is  given  in  Chapter  V. 
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III.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  experimental  studies  were  carried  out  in  two  phases. 
Phase  I  consisted  of  studies  employing  a  homogeneous  mixture 
combustion  of  methane  and  acetylene  while  in  Phase  II  a  non- 
homogeneous  mixture  combustion  was  used.  The  various  tests 
conducted  are  summarized  in  Table  3.  1.  The  tests  in  Phase  I 
were  conducted  in  five  parts,  of  which  (a)  and  (b)  were  explana¬ 
tory  tests.  The  tests  in  Phase  II  were  conducted  in  three  parts. 
These  are  described  in  more  detail  later;  however,  as  noted, 
these  tests  were  conducted  to  determine  the  effect  of  the  com¬ 
bustion  of  various  mixtures  and  conditions  of  combustion  on 
combustion  time  and  the  trajectory  of  the  flame. 

The  combustion  time  is  defined  in  the  present  work  as  the 
interval  between  the  time  of  ignition  and  the  time  the  combustion 
pressure  reached  its  maximum  value.  The  combustion  time  was 
adopted  as  one  of  the  important  factors  in  the  present  study  for  the 
following  reason.  From  the  definition  of  heat  of  combustion,  it 
can  be  shown  that  the  higher  the  time  rate  of  the  combustion 
pressure  is,  the  higher  the  time  rate  of  energy  emission  from  the 
combustion  becomes .  Therefore,  for  a  constant  maximum  com¬ 
bustion  pressure,  the  time  rate  of  the  combustion  pressure 
becomes  almost  inversely  proportional  to  the  combustion  time  as 
defined  above.  During  the  initial  stage  of  the  experimental  investi¬ 
gation,  the  maximum  combustion  pressure  was  found  to  be  constant 
as  long  as  the  amount  of  flammable  gas  mixture  charged  into  the 
combustion  chamber  was  fixed.  Consequently  it  was  decided 
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to  use  combustion  time  as  a  parameter  indicating  the  time  rate  of 
energy  emission  which  undoubtedly  is  one  of  the  most  important 
factors  in  the  study  of  combustion. 

3.  1  Apparatus  and  Experimental  Procedure  for  Phase  I. 

3.1.1  Mechanical  Apparatus 

In  Phase  I,  an  investigation  was  made  of  the  effect  of 
combustion  chamber  rotative  speed  and  location  of  the  point  of 
ignition  on  combustion  time  and  maximum  pressure  and  of  the 
flame  trajectories.  In  order  to  determine  accurately  the  flame 
trajectory  by  experimental  methods,  the  equipment  was  designed 
so  that  all  conditions  in  the  combustion  chamber  would  be  as  true 
as  possible  to  the  assumptions  made  for  the  analytical  solution 
of  the  flame  trajectory,  namely,  solid  body  rotation  and  negligible 
wall  effect.  Figure  3.  1  is  a  schematic  of  the  setup  used  for  two 
parts  of  the  investigation  of  the  homogeneous  mixture  combustion. 
The  first  part  concerns  the  effect  of  the  combustion  chamber 
rotative  speed  and  location  of  the  point  of  ignition  on  both  com¬ 
bustion  time  and  pressure,  and  the  second  part  concerns  the  flame 
trajectory. 

Combustion  Chamber: 

Since  only  a  small  amount  of  methane -air  or  acetylene -air 
mixture  was  injected  into  the  combustion  chamber  at  room  tempe¬ 
rature  and  atmospheric  pressure,  and  since  combustion  was  effected 
only  intermittently,  the  strength  of  the  combustion  chamber  wa»  no 
serious  problem.  During  the  conduct  of  the  test,  the  maximum 


FIG.  3.1.  EXPERIMENTAL  SETUP  FOR  THE  HOMOGENEOUS  MIXTURE  COMBUSTION 
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chamber  pressure  encountered  was  approximately  140  pounds  per 
square  inch  and  the  maximum  chamber  temperature  approximately 
230  degrees  Fahrenheit. 

Figure  3.  2  shows  a  cross-section  of  the  combustion  chamber 
used  for  the  homogeneous  mixture  combustion  studies.  The 
combustion  chamber  was  made  of  ordinary  steel,  cylindrical  in 
shape,  2.  5  inches  inside  diameter  and  0.  5  inches  long.  One  end 
of  the  combustion  chamber  was  made  of  Luc  it  e  so  that  the  flame 
could  be  photographed  with  a  Fastax  camera.  Aniline,  which  was 
added  to  the  fuel  to  increase  the  brightness  of  the  flame,  reacts 
chemically  with  Lucite;  therefore,  a  piece  of  glass  was  placed  over 
the  Lucite  O'"  inside  of  the  chamber.  In  order  to  prevent 
leakage  around  the  glass  window,  a  rubber  gasket  was  inserted 
between  the  combustion  chamber  and  the  glass  window;  and  the 
joint  around  the  glass  and  Lucite  was  sealed  with  putty.  The 
other  end  of  the  chamber  was  steel  and  had  a  smooth  counterbored 
hole  for  the  insertion  of  a  pressure  pickup  mounted  in  a  fixed 
holder  as  shown  in  Fig.  3.  2. 

Figure  3.2  is  a  photograph  of  an  exploded  view  of  the  com¬ 
bustion  chamber.  Spark  plugs  and  valves  were  symmetrically 
inserted  to  provide  for  good  balance  as  the  chamber  was  rotated 
at  speeds  up  to  3600  rpm.  Only  one  of  the  spark  plugs  was  used, 
however,  while  the  others  were  inserted  for  balancing  only.  Since 
the  combustion  chamber  was  small,  the  combustion  pressure  was 
essentially  uniform  throughout  the  entire  chamber;  therefore,  the 
indicated  pressure  at  the  center  of  the  combustion  chamber  where 


Contact  Brush 


FIG.  3.  2  CROSS-SECTION  OF  THE  COMBUSTION  CHAMBER  (FULL  SCALE) 
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the  pressure  pickup  was  located  represented  the  chamber  pressure. 

The  combustion  chamber  was  threaded  unto  a  hollow  steel 
cylinder  as  shown  in  Fig.  3.  2.  The  hollow  cylinder  was  supported 
by  a  journal  bearing  that  was  designed  for  a  main  bearing  of  a 
G.M.  C.  automobile  crankshaft.  Since  the  operating  time  for  the 
combustion  chamber  was  only  a  few  minutes,  the  journal  bearing 
was  lubricated  by  use  of  a  gravity  lubricator  with  SAE  30  engine 
oil.  Thus  the  hollow  cylinder  provided  a  rotational  mounting  for 
the  combustion  chamber.  A  V-groove  on  the  hollow  cylinder  was 
belted  to  a  variable  speed  drive  driven  by  an  electric  motor.  The 
cylinder  with  the  combustion  chamber  attached  could  thus  be  driven 
at  speeds  from  0  to  3600  rpm. 

Spark  Plug: 

The  spark  plugs  used  were  AC  type  M-8  (Fig.  3.  3).  The 
active  spark  plug  was  connected  to  a  copper  ring  supported  by  a 
Bakelite  ring  mounted  on  the  outside  of  the  combustion  chamber. 
This  ring  provided  an  electrical  connection  from  the  induction  coil 
(Fig.  3.  1  and  3.  2)  to  the  rotating  chamber.  The  electrode  of  the 
active  spark  plug  was  modified  with  copper  wires  of  different 
lengths  to  provide  for  ignition  at  three  points  in  the  cylindrical 
chamber --at  the  periphery,  point  of  half  radius  from  the  center, 
and  at  the  center. 

Pressure  Pickup: 

The  pressure  was  measured  using  a  piezoelectric  pickup 
with  a  range  of  0  to  1000  psi,  electrically  connected  to  a  Model 
533  Tektronix  oscilloscope  through  a  calibrator.  The  diaphragm 


FIG.  3«3 .  COMBUSTION  CHAMBER 
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end  of  the  pressure  pickup  was  threaded  into  an  adaptor  which  was 
supported  by  the  holder  as  shown  in  Fig.  3.  2.  Thus  the  pressure 
pickup  was  stationary,  and  as  noted,  the  combustion  chamber  was 
rotatable.  To  prevent  leakage  around  the  contact  surface  between 
the  pressure  pickup  adaptor  and  the  combustion  chamber,  a  Teflon 
cap  was  attached  to  the  adaptor  with  epoxy  cement.  This  provided 
an  essentially  leak  free  bearing  surface  between  the  pressure 
pickup  and  the  rotative  combustion  chamber,  with  force  applied 
from  the  pressure  pickup  holder. 

Driving  Motor: 

An  A-C  induction  motor  with  a  variable  speed  drive  was 
connected  by  a  single  V-belt  to  the  combustion  chamber  as  shown 
in  Fig.  3.  1.  This  motor  was  rated  1/2  hp,  1760  rpm,  222  volts, 
single  phase. 

Premixing  Tank: 

In  parts  (a)  and  (c)  of  Phase  I,  methane  was  used  to 
determine  the  effect  of  rotative  combustion  chamber  speed  and 
location  of  the  point  of  ignition  on  combustion  time  and  maximum 
combustion  pressure.  (As  noted  later,  acetylene  was  used  to 
determine  the  flame  trajectory.) 

Before  the  methane -air  mixture  was  charged  to  the  com¬ 
bustion  chamber,  the  methane  and  air  were  introduced  to  a  pre- 
mixing  tank  where  they  could  become  mixed  homogeneously.  The 
mixture  ratio  of  methane  to  air  could  also  be  controlled  in  this 
tank.  The  premixing  tank  was  made  of  plastic  with  1/2  Inch 
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scale  markings  on  the  surface.  The  diameter  of  the  tank  was  6 
inches  and  the  height  was  7.  5  inches.  Thus,  the  volume  of  gas 
charged  to  the  combustion  chamber  was  determined  by  measuring 
on  the  scale  markings  the  displaced  water  in  the  premixing  tank. 

The  accuracy  of  controlling  the  mixture  ratio  in  this  manner  was 
adequate  for  the  purpose  of  these  tests. 

In  the  determination  of  the  flame  trajectory  with  acetylene, 
the  weight  of  the  actylene-air  mixture  had  to  be  more  precisely 
measured,  since  the  flame  propagation  speed  of  the  acetylene  - 
air  mixture  is  very  sensitive  to  variations  of  the  mixture  ratio  (14). 
Therefore,  the  premixing  tank  was  changed  to  a  steel  tank  (of  21 
cubic  inch  volume)  which  was  connected  to  a  77  inch  manometer 
by  a  1/4  inch  copper  tube.  More  precise  measurements  of  mixture 
ratios  were  possible  with  the  manometer.  This  steel  premixing 
tank  was  also  used  for  the  nonhomogeneous  mixture  combustion 
tests . 

3.  1.  2  Photographic  Equipment 

As  shown  in  Fig.  3.  1  the  oscilloscope  was  connected  to  the 
piezocalibrator  as  an  indicator  to  show  the  pressure  variation 
during  combustion.  A  Polaroid  camera,  using  Polaroid  Type  42 
film,  was  attached  to  the  oscilloscope  in  order  to  make  a  permanent 
record  of  the  measurements  of  pressure  variation.  The  oscillo¬ 
scope  was  intensity  modulated  from  a  time  mark  generator  to  give 
an  accurate  time  reference  on  the  photographs. 

High  speed  photographs  of  the  flame  during  combustion 
were  made  through  the  Lucite  and  glass  window  on  the  end  of  the 
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combustion  chamber  in  order  to  determine  experimentally  the 
flame  trajectory  and  in  order  to  observe  flame  spread.  These 
photographs  were  taken  with  a  16mm  Fastax  camera,  with  its 
lens  located  28  inches  from  the  center  of  the  combustion  chamber 
and  the  aperture  fully  open.  As  noted  in  Fig.  3.  1  and  3.4,  the 
photographs  were  taken  through  the  center  of  the  hollow  steel 
cylinder  supporting  the  combustion  chamber.  Since  the  combustion 
time  was  of  the  order  of  10  milliseconds,  only  about  25  feet  of 
Kodak  Tri-x  film  was  used  for  each  record.  This  camera  produced 
a  1/120  second  time  mark  on  the  film  so  that  an  accurate  time 
scale  could  be  obtained  for  the  flame  trajectory  photographs, 
independent  of  the  small  variations  in  the  camera  speed. 

3.  1.  3  Electrical  System: 

Since  this  investigation  was  to  determine  the  various  effects 
during  combustion  time,  and  since  the  physical  phenomena  occur 
during  a  period  of  about  5  to  10  milliseconds,  an  electric  recording 
system,  synchronized  with  the  normal  electrical  ignition  system, 
was  necessary. 

The  ignition  system  and  the  recording  system  are  shown  in 
schematic  form  in  Fig.  3.  5.  The  ignition  system  consisted  of  a  6 
volt  battery,  switches  SI  and  S2,  solenoid  relay,  an  induction  coil, 
and  a  spark  plug.  The  recording  system  consisted  of  the  piezo- 
calibrator,  Tektronix  type  181  time  mark  generator,  audio 
frequency  amplifier,  and  Tektronix  oscilloscope  type  533.  The 
time  mark  generator  produced  1  millisecond  time  marks  by 
intensity  modulation  of  the  oscilloscope. 


FIG.  3.  5.  WIRING  DIAGRAM  FOR  EXPERIMENTAL  SETUP 
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The  piezo  calibrator  provided  a  calibration  voltage  so  that 
the  vertical  deflection  of  the  oscilloscope  could  be  accurately 
calibrated  in  psi  to  the  signal  from  the  piezoelectric  pickup. 

The  type  533  Tektronix  oscilloscope  has  a  rise  time  of 
0.025  microseconds,  which  is  more  than  adequate  to  record 
accurately  the  pressure  variations  which  occurred  during  the 
10  millisecond  combustion  time. 

3.1.4  Experimental  Procedure: 

In  parts  (a)  and  (c)  of  Phase  I,  methane  gas  was  used  to 
determine  the  effect  of  the  combustion  chamber  rotative  speed 
and  location  of  the  point  of  ignition  on  the  combustion  time  and 
chamber  pressure.  Methane  gas  was  introduced  into  the  pre - 
mixing  tanks  and  measured  by  the  volume  of  water  displaced. 

Air  was  then  introduced  into  the  premixing  tank  until  the  methane 
to  air  ratio  was  1:10  in  volume.  In  order  to  insure  a  homogeneous 
mixture  of  air  and  methane,  the  mixture  was  kept  in  the  premixing 
tank  for  a  few  minutes. 

The  combustion  chamber  was  scavenged  by  the  methane-air 
mixture,  and  then  a  charge  of  the  methane -air  mixture  was  intro¬ 
duced  into  the  combustion  chamber.  The  volume  of  the  mixture 
charged  was  three  times  that  of  the  chamber  volume  at  atmospheric 
pressure.  The  volume  of  the  charge  was  measured  by  the  scale 
marked  on  the  surface  of  the  premixing  tank. 

During  the  test  with  methane,  it  was  discovered  that  the 
color  of  the  methane  flame  was  unsuited  to  the  photographic 
recording  of  the  flame  trajectory.  Also,  the  ratio  of  the  density 
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of  the  air  to  that  of  the  burning  methane -air  mixture  flame  is 
unknown,  but  this  ratio  was  needed  to  calculate  the  flame  trajectory 
as  shown  by  Equation  2.  9  and  2.  10.  The  acetylene-air  mixture 
with  a  few  drops  of  aniline  added  proved  to  give  quite  satisfactory 
flame  photographs;  and  the  density  ratio  of  an  air  to  acetylene -air 
mixture  flame  is  known  to  be  approximately  3.0.  Therefore,  in 
parts  (b)  and  (d)  of  Phase  I,  acetylene  was  substituted  in  order 
to  obtain  all  needed  data  in  connection  with  the  flame  trajectories. 

As  noted  in  Section  3.  1.  1,  the  steel  mixing  tank  with  the 
77  inch  manometer  for  more  precise  measurements  was  used  for  the 
acetylene.  This  premixing  tank,  with  air  at  atmospheric  pressure, 
was  then  charged  with  acetylene  gas  to  a  pressure  of  19  inches  of 
mercury  above  atmospheric  pressure.  This  resulted  in  an  acetylene 
to  air  mixture  ratio  of  IT  1.74  in  weight.  This  mixture  was  allowed 
to  remain  in  the  tank  until  the  air  remaining  in  the  mixing  tank 
could  be  assumed  to  have  mixed  homogeneously  with  the  acetylene. 
After  scavenging  the  combustion  chamber  with  air,  the  acetylene - 
air  mixture  was  then  charged  into  the  combustion  chamber  from 
the  premixing  tank  through  the  inlet  valve  (1/8  inch  pipe  needle 

valve)  until  the  pressure  in  the  premixing  tank  decreased  by  1.  5 

"4 

inches  of  mercury  (0.  316  x  10  lb). 

During  all  tests  in  Phase  I,  before  ignition  was  effected, 
the  A-C  motor  with  the  variable  speed  drive  was  switched  on  and 
the  speed  of  the  combustion  chamber  was  held  constant  for  one 
or  two  minutes  in  order  to  ensure  that  the  gas  mixture  within  the 
combustion  chamber  attained  a  constant  speed  of  solid  body 
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rotation  on  one  hand  (7)  and  that  the  fuel -air  mixture  became 
mixed  homogeneously  on  the  other  hand.  Next,  switch  SI  was 
turned  on  for  a  few  seconds  to  excite  the  induction  coil,  the  shutter 
of  the  Polaroid  camera  was  opened,  the  Fastax  camera  and  switch 
S2  were  turned  on.  Upon  completion  of  the  above  procedure,  at 
the  end  of  combustion,  the  shutter  of  the  Polaroid  camera  was 
closed.  Figure  3.6  shows  the  flow  diagram  of  this  procedure. 
When  switch  SI  is  turned  on,  the  induction  coil  is  energized. 

When  switch  S2  is  turned  on,  the  solenoid  relay  is  opened  and  the 
secondary  of  the  induction  coil  is  brought  to  a  high  voltage  which 
causes  the  spark  in  the  combustion  chamber.  This  high  ignition 
voltage  also  initiates  the  horizontal  sweep  on  the  oscilloscope 
through  the  triggering  input  and  produces  a  spike  on  the  oscillo¬ 
scope  trace.  The  Polaroid  camera  photographed  this  spike, 
which  established  the  beginning  of  combustion.  Since  the  trace 
on  the  oscilloscope  was  also  marked  with  1  millisecond  time 
marks,  the  relationship  between  combustion  pressure  and  time 
and  total  combustion  time  was  easily  found. 

Since  the  combustion  chamber  was  small,  combustion 
pressure  in  the  entire  chamber  was  assumed  to  be  uniform. 
Therefore  the  curve  of  the  relationship  between  combustion 
pressure  and  time,  which  was  recorded  by  the  Polaroid  camera, 
may  be  assumed  to  represent  the  relationship  between  the  com¬ 
bustion  pressure  and  time  for  the  entire  combustion  chamber. 

Following  the  procedure  described  above,  the  effect  of 
the  rotative  chamber  speed  and  location  of  the  point  of  ignition 
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FIG.  3.  6.  FLOW  DIAGRAM  OF  EXPERIMENTAL  PROCEDURE 
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on  combustion  time  and  maximum  combustion  pressure  were  made 
at  rotative  combustion  chamber  speeds  of  0,  500,  900,  1350,  1800, 
2700  and  3600  rpm.  Throughout  the  experiments,  the  combustion 
chamber  was  charged  with  a  methane -air  mixture  of  three  times 
the  volume  of  the  chamber,  and  a  constant  methane -air  mixture 
ratio  of  1:10  was  maintained. 

For  part  (d)  of  Phase  I,  in  the  experimental  determination 
of  flame  trajectories,  observations  were  made  at  the  rotative  com¬ 
bustion  chamber  speed  of  900  and  1800  rpm  with  ignition  taking 
place  at  a  point  of  half  radius  from  the  center.  An  acetylene -air 

mixture  with  a  mixture  ratio  of  1:1.  74  and  an  amount  of  0.  316  x 
-4 

10  lb  was  used  throughout  the  tests, 

3.1.5  Flame  Photography: 

The  flame  color  of  the  methane  gas  mixture  is  blue  at 
any  mixture  ratio,  and  is  rather  difficult  to  photograph.  Therefore 
as  noted  before,  an  acetylene  mixture  was  used  for  the  experimental 
determination  of  the  flame  trajectory  because  it  gave  a  red  color 
flame  at  a  rich  mixture  ratio,  thus  allowing  sufficient  light  for  film 
exposure.  Carbon  particles  which  were  deposited  on  the  glass 
window  as  a  result  of  combustion  were  removed  after  each  experi¬ 
ment.  Even  with  acetylene,  a  lean  mixture  ratio  resulted  in  a 
red/blue  flame  which  was  difficult  to  photograph.  To  obtain 
better  flame  photographs  in  this  case,  it  was  necessary  to  intensify 
the  color  of  the  flame.  For  this  purpose,  a  catalyst,  aniline 
(C^Hj-C^)*  was  used.  Aniline  contains  a  CN  radical  which,  when 
burning,  emits  an  orange  color  which  is  strong  enough  to  record 
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on  Kodak  Tri-x  film.  In  order  to  achieve  a  homogeneous  mixture, 
a  few  drops  of  aniline  were  added  to  the  acetylene -air  mixture  in 
the  steel  premixing  tank,  thus  allowing  evaporated  aniline  to  mix 
with  the  fuel  mixture.  Other  catalysts  which  could  have  been  used 
are  metyl  amine  (CH^NH^)*  metyl  nitrate  (CH^NO^),  or  nitrobenzen 
(C^HgNO^).  Of  these,  only  nitrobenzen  was  tried,  but  it  was  found 
that  aniline  gave  better  results.  Since  aniline  is  very  toxic  it  was 
necessary  to  avoid  inhaling  the  fumes. 

3.  2  Apparatus  and  Experimental  Procedure  for  Phase  II. 

In  Phase  II  of  the  investigation  the  premixing  steel  tank, 
the  ignition  system,  and  recording  system  were  the  same  as  those 
used  in  Phase  I.  In  this  section  the  modified  combustion  chamber, 
fuel  charging  system,  and  the  necessary  photographic  technique 
will  be  described.  An  attempt  has  been  made  to  develop  an 
empirical  method  of  predicting  the  effect  of  the  combustion  chamber 
rotative  speed  and  location  of  the  point  of  ignition  on  the  combustion 
time  for  the  nonhomogene ous  mixture  combustion. 

3.  2.  1  Modified  Combustion  Chamber  and  Charging  System 

In  order  to  produce  a  nonhomogeneous  mixture  in  the  com¬ 
bustion  chamber,  the  chamber  was  modified  to  simulate  the 
conditions  of  a  diesel  engine.  The  combustion  chamber  was 
modified  as  shown  in  Fig.  3.  7.  The  modification  was  made  oply 
on  the  pressure  pickup  side  to  provide  a  means  of  charging  the 
combustion  chamber  while  it  was  being  rotated.  A  small  copper 
tube  1/16  inch  in  diameter  was  connected  directly  to  the  diaphragm 
side  of  the  pressure  pickup  through  the  adaptor  utilizing  the 
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FIG.  3.  8.  EXPERIMENTAL  SETUP  FOR  NONHOMOGENEOUS  MIXTURE  COMBUSTION. 
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clearance  around  the  tube  as  an  inlet  passage  into  charging  ports 
drilled  in  the  chamber  cover.  Charging  ports  were  made  at  the 
periphery,  half  radius  from  the  center,  and  the  center  of  the 
combustion  chamber.  The  desired  charging  port  was  left  open 
while  the  other  two  were  closed  during  an  experiment.  The 
location  of  the  point  of  ignition  was  also  changed  to  one  of  the 
three  positions  (periphery,  half  radius  from  the  center,  or  center) 
for  each  of  the  charging  locations.  This  modification  enabled  the 
observation  of  flame  spread  to  be  made  with  any  combination  of 
rotative  combustion  chamber  speed,  location  of  charging  port, 
and  location  of  the  point  of  ignition. 

In  the  case  of  the  flame  trajectory  measurements  it  was 
necessary  to  synchronize  the  charging  of  the  fuel  to  the  com¬ 
bustion  chamber  with  the  ignition  and  recording  events.  This 
need  for  synchronizing  is  due  to  the  fact  that  it  is  desirable,  as 
explained  in  5.  2,  to  ignite  the  flammable  mixture  while  it  is 
still  contained  within  the  restricted  region  around  the  charging 
port,  instead  of  having  it  spread  throughout  the  chamber.  A 
solenoid  valve  was  inserted  in  the  line  between  the  steel  premixing 
tank  and  modified  combustion  chamber  to  accomplish  this 
synchronization. 

3.  2.  2  Experimental  Procedure; 

In  order  to  find  the  effects  of  the  combustion  chamber 
rotative  speed,  location  of  the  point  of  ignition,  and  charging  port 
position  on  combustion  time  and  on  the  flame  trajectory  for  the 
nonhomogeneous  mixture  combustion,  parts  (f) ,  (g),  and(h),  the 


modified  combustion  chamber  and  charging  system  were  used. 

An  overrich  acetylene -air  mixture  was  used  in  this  part  of  the 
nonhomogeneous  mixture  combustion  tests. 

The  acetylene  gas  was  charged  into  the  steel  premixing 

tank  until  the  pressure  was  10  inches  of  mercury  above  atmospheric. 

Then  air  was  charged  to  the  tank  until  the  absolute  pressure  was 

43.92  inches  of  mercury  above  atmospheric  (the  acetylene -air 

mixture  ratio  was  1:4.85  in  weight).  When  the  combustion 

chamber  reached  the  desired  speed,  the  acetylene -air  mixture 

was  charged  into  the  chamber  through  the  desired  charging  port 

(refer  to  Fig.  3.  7  and  3.  8)  until  the  pressure  in  the  premixing  tank 
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decreased  to  1  inch  of  mercury  (0. 3  x  10  lb  in  weight).  In  the 
investigation  of  the  variation  in  combustion  time,  parts  (f)  and  (g), 
it  is  necessary  to  have  the  mixture  ignited  soon  after  the  charging. 
However,  there  is  no  need  to  achieve  such  a  critical  synchroni¬ 
zation  between  the  charging  of  the  mixture  and  the  ignition  as 
necessary  in  the  study  of  flame  trajectories,  part  (h) ,  since  the 
condition  of  nonhomogeneity  of  the  mixture  persists  in  the  chamber 
for  a  reasonable  period  of  time.  Therefore  no  attempt  was  made 
to  synchronize  the  charging  with  ignition  and  recording  for  the 
measurements  in  parts  (f)  and  (g).  The  charging  valve  used  was 
a  1/4  inch  pipe  needle  valve,  hand  operated.  The  firing  and 
recording  procedures  were  the  same  as  outlined  in  Section  3.  lf  3. 

In  the  experimental  determination  of  the  flame  trajectory 
it  was  necessary  to  control  accurately  the  quantity  of  the  fuel 
charged  to  the  combustion  chamber,  and  to  synchronize  the  ignition 
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with  fuel  charging  and  the  recording  events.  The  synchronization 
was  accomplished  by  the  use  of  the  solenoid  valve.  However, 
since  the  pressure  at  the  charging  port  changed  with  combustion 
chamber  speed,  the  amount  of  gas  mixture  entering  the  chamber 
when  the  solenoid  valve  was  opened  varied  depending  on  the  com¬ 
bustion  chamber  rotative  speed  and  thus  made  it  rather  difficult  to 
control  the  amount  of  mixed  fuel  that  was  charged. 

It  was  experimentally  determined  that,  when  acetylene  was 
charged  into  the  steel  premixing  tank  which  was  originally  filled 
with  air  to  a  pressure  of  3  inches  of  mercury,  an  acetylene -air 
mixture  with  a  mixture  ratio  of  1:11  resulted  in  the  tank.  Further¬ 
more,  when  the  solenoid  valve  was  opened  to  charge  the  mixture 
into  the  combustion  chamber  which  was  rotating  at  900  rpm, 

the  pressure  in  the  tank  was  allowed  to  decrease  by  0.  3  inch  of 
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mercury,  indicating  that  0.09  x  10  pounds  of  mixture  had  been 
charged  into  the  chamber. 

In  view  of  such  a  lean  mixture  in  the  combustion  chamber, 
it  was  necessary  to  intensify  the  brightness  of  the  flame  for  flame 
photography  by  adding  a  few  drops  of  aniline  to  the  steel  premixing 
tank  as  noted  previously  in  Section  3.  1.4.  After  the  combustion 
chamber  obtained  a  speed  of  900  rpm,  switch  SI  (refer  to  Figures 
3.  5  and  3.  6)  was  turned  on  and  the  shutter  of  the  Polaroid  camera 
was  opened.  The  solenoid  valve,  the  Fastax  camera,  and  switch 
S2  were  turned  on  simultaneously  (refer  to  Figures  3.  5  and  3.6). 
After  combustion  was  completed,  the  shutter  of  the  Polaroid 


camera  was  closed. 
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As  in  Phase  I  for  the  homogeneous  case  and  using  the  same 
procedure,  photographic  recordings  were  taken  of  the  combustion 
pressure  as  a  function  of  time  on  the  oscilloscope  with  the  Polaroid 
camera  and  of  the  flame  itself  with  the  Fastax  camera.  The 
results  of  the  experimental  work  are  described  in  Chapter  IV. 
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IV.  RESULTS  OF  EXPERIMENT  AND 

— THEORETICAL  CALCULATTOTT — 

4.  1  Results  of  Preliminary  Tests 

Two  preliminary  tests  were  carried  out  for  the  purpose  of 
observing  the  general  behavior  of  the  combustion  pressure  variation 
and  the  development  of  flame  during  combustions. 

For  the  first  test  part  (a)  of  Phase  I,  a  flammable  mixture 
of  methane  gas  and  air  with  a  relatively  slow  burning  velocity  was 
used.  Combustion  pressure  variations  were  obtained  and  compared 
for  chamber  rotative  speeds  of  1800  and  2700  rpm  when  the  mixture 
was  ignited  at  the  periphery  of  the  chamber.  These  results  are 
summarized  in  Fig.  4.  1,  from  which  it  can  be  seen  that  the  maxi¬ 
mum  combustion  pressure  remained  reasonably  constant  in  spite 
of  the  fact  that  the  combustion  time  varied.  Variation  in  com¬ 
bustion  time  is  very  large  at  2,  700  rpm  as  indicated  in  Fig.  4.  1. 

This  is  believed  to  have  been  caused  by  vibration  of  the  combustion 
chamber.  During  the  experiments  the  chamber  was  frequently 
observed  to  resonate  around  2,  700  rpm. 

For  the  second  test,  part  (b)  of  Phase  I,  a  Fastax  high 
speed  camera  was  used  in  order  to  photograph  the  flame  develop¬ 
ment  under  the  influence  of  centrifugal  and  drag  forces.  High  speed 
pictures  of  the  flame  were  taken  using  a  rich  acetylene -air  mixture, 
which  burns  with  high  flame  luminescense.  It  was  observed  that  at 
high  chamber  rotative  speeds  the  flame,  when  ignited  at  the  pe¬ 
riphery  of  the  chamber,  moved  radially  inward  describing  a  spiral 
curve.  Afterward,  the  flame  moved  outward  from  the  center  until 


the  entire  chamber  was  filled  with  the  flame. 
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FIG.  4.1.  CHAMBER  PRESSURE  VS.  COMBUSTION  TIME  FOR 
VARIOUS  ROTATIVE  CHAMBER  SPEED  FOR  THE 
HOMOGENEOUS  MIXTURE  COMBUSTION. 


49 


It  is  noteworthy  that  a  pair  of  appreciably  strong  swirls  of 
opposite  sense  resulted  as  the  flame  disappeared,  as  shown  in 
Fig.  4.  2  and  Fig.  4.  3.  In  the  static  chamber  the  flame  propagation 
is  observed  to  have  the  shape  shown  in  Fig.  4.  4. 

4.  2  Results  of  the  Effect  of  Location  of  the  Point  of  Ignition  and 

Chamber  Rotative  Speed  on  Combustion  Time  for  Homogeneous 
Mixture  Combustion. 

The  combustion  times  were  compared  for  combustion 
chamber  rotative  speeds  of  0,  500,  900,  1300,  1800,  2700,  and 
3600  rpm  and  at  three  locations  of  the  point  of  ignition:  at  the 
periphery,  half  radius  from  the  center,  and  the  center.  Through¬ 
out  the  tests,  the  combustion  chamber  was  charged  with  methane - 
air  mixture  as  specified  in  Part  (c)  of  Table  3.  1.  The  results 
are  summarized  in  Fig. 4.  5. 

In  the  case  of  ignition  at  the  periphery,  the  influence  of  the 
combustion  chamber  rotative  speed  on  the  combustion  time  is  most 
distinct.  The  combustion  time  at  high  speeds  was  approximately 
40%  lower  than  the  combustion  time  for  the  stationary  chamber. 

In  the  case  of  ignition  at  the  center,  on  the  other  hand, 
the  combustion  time  was  less  influenced  by  the  combustion  chamber 
rotative  speed.  The  combustion  time  was  almost  independent  of 
the  chamber  speed. 

4.  3  Results  of  the  Effect  of  Location  of  the  Point  of  Ignition  and 
Combustion  Chamber  Rotative  Speed  on  Combustion  Time  for 
the  Nonhomogeneous  Mixture  Combustion. 

The  effect  of  changes  in  the  location  of  the  point  of  ignition 


FIG.  4.  2  FLAME  PATTERN  IN  HOMOGENEOUS  MIXTURE 
COMBUSTION  AT  2400  RPM. 

4500  FRAMES  PER  SECOND 
MIXTURE  RATIO  (ACETYLENE  TO  AIR)  1:7 
(CONTINUOUS  SINGLE  FILM  STRIP  DISPLAYED  IN 
SIX  SECTIONS,  A  THROUGH  F. ) 


FIG.  4.  3.  FLAME  PATTERN  IN  HOMOGENEOUS  MIXTURE 
COMBUSTION  AT  2700  RPM. 

4500  FRAMES  PER  SECOND. 

MIXTURE  RATIO  (ACETYLENE  TO  AIR)  1:7 
(CONTINUOUS  SINGLE  FILM  STRIP  DISPLAYED 
IN  FOUR  SECTIONS,  A  THROUGH  D. ) 


FIG.  4.  4.  FLAME  PATTERN  IN  HOMOGENEOUS  MIXTURE 
COMBUSTION  AT  STATIONARY. 

4500  FRAMES  PER  SECOND. 

MIXTURE  RATIO  (ACETYLENE  TO  AIR)  1:7 
(CONTINUOUS  SINGLE  FILM  STRIP  DISPLAYED 
IN  FIVE  SECTIONS,  A  THROUGH  E. ) 
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FIG.  4.  5.  COMBUSTION  TIME  VS.  ROTATIVE 
CHAMBER  SPEED  FOR  THE  THREE 
IGNITION  POSITIONS.  (HOMOGENEOUS 
CASE). 
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and  of  the  combustion  chamber  rotative  speed  on  combustion  time 
for  the  nonhomo  gene  ous  mixture  combustion  was  investigated,  using 
the  apparatus  shown  in  Fig.  3.  7  and  Fig.  3.  8,  for  three  locations 
of  the  point  of  ignition  (periphery,  half  radius  from  the  center, 
and  center)  and  for  two  locations  of  the  charging  port,  one  at  the 
periphery  and  the  other  at  half  radius  from  the  center.  Through¬ 
out  these  tests,  an  overrich  acetylene -air  mixture,  as  specified 
in  Part  (f)  of  Table  3.  1,  was  admitted  through  the  charging  port. 

The  results  of  these  tests  are  shown  in  Fig.  4.  6  and  4.  7.  Each 
point  given  in  Fig.  4.  6  and  Fig.  4.7  represents  the  arithmetic 
mean  value  of  four  tests  carried  at  each  chamber  rotative  speed. 

The  experiment  was  first  carried  out  using  the  charging 
port  at  the  periphery.  In  the  case  of  ignition  at  the  periphery, 
the  influence  of  combustion  chamber  rotative  speed  on  the  com¬ 
bustion  time  is  evident  in  Fig.  4.  6.  The  combustion  time  at  high 
speeds  was  approximately  30%  lower  than  the  combustion  time 
when  the  chamber  was  stationary.  In  the  case  of  ignition  at  the 
center  (refer  to  the  description  in  3.  1.1),  the  combustion  time 
was  influenced  little  by  combustion  chamber  rotative  speeds;  the 
combustion  time  at  high  speeds  was  approximately  only  5%  lower 
than  that  when  the  chamber  was  stationary. 

In  the  second  experiment  of  Part  (f) ,  the  acetylene -air 
mixture  was  charged  at  a  point  of  half  radius  from  the  center 
(refer  to  Fig.  3.  7).  and  the  three  locations  of  the  point  of  ignition 
(periphery,  half  radius  from  the  center,  and  center)  were  used. 

As  shown  in  Fig.  4.  7,  the  combustion  time  increased  slightly  as 
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the  chamber  speed  increased  for  all  three  locations  of  the  point 
of  ignition.  This  increase  in  combustion  time  was  about  5%  of 
that  when  the  chamber  was  "stationary.  When  ignition  was  effected 
at  each  typical  position  (periphery,  half  radius  from  the  center, 
and  center),  a  pressure  spike  was  recorded  at  the  beginning  of 
combustion  on  the  photographs  which  were  taken  by  the  Polaroid 
camera.  The  reason  for  this  photography  will  be  discussed  in 
Chapter  V. 

4.4  Comparison  of  the  Effect  of  Combustion  Chamber  Rotative 
Speed  on  Combustion  Time  for  Homogeneous  and  Non- 
homogeneous  Mixture  Combustion. 

For  the  purpose  of  investigating  the  difference  in  the 
behavior  of  combustion  in  the  homogeneous  and  nonhomogeneous 
mixture  combustion  cases,  the  effect  of  combustion  chamber 
rotative  speed  on  combustion  time  was  studied  for  the  two  cases 
using  overrich  acetylene  to  air  mixtures  as  specified  in  Part  (e)  and 
(g)  of  Table  3.  1.  The  results  shown  in  Fig.  4.  8  indicate  that  the 
difference  in  the  combustion  time  between  the  homogeneous  and 
the  nonhomogeneous  mixture  combustion  is  large. 

4.  5  Experimental  Determination  of  Flame  Trajectory. 

In  order  to  determine  experimentally  the  flame  tra¬ 
jectories  and  measure  the  flame  radii,  photographs  were  taken 
with  the  Fastax  camera  at  900  and  1800  rpm  for  the  case  of 
homogeneous  mixture  combustion  with  the  ignition  effected  at 
half  radius  from  the  center  of  the  combustion  chamber.  These 
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Rotative  Chamber  Speed  (rpm) 

FIG.  4.8.  COMPARISON  OF  COMBUSTION  TIME  FOR  THE 
HOMOGENEOUS  AND  NONHOMOGENEOUS  OVER¬ 
RICH  MIXTURE. 


pictures,  which  are  shown  in  Fig.  4.9  and  4.  10,  were  taken 
with  the  acetylene -air  mixture  as  specified  in  Part  (d)  of  Table 
3.  1. 

For  analysis  these  photographs  were  enlarged  to  1.  16 
times  actual  size.  The  Fastax  camera  enabled  the  sparking 
point,  that  is  the  beginning  point  of  combustion,  to  be  photo¬ 
graphed.  From  these  photographs,  which  distinctly  showed  the 
outside  contour  of  the  combustion  chamber,  the  center  of  the 
combustion  chamber  was  also  obtained. 

It  was  rather  difficult  to  determine  the  center  of  the 
flame  because  the  flame  did  not  form  a  perfect  circle.  It  was 
therefore  necessary  to  approximate  such  a  circle  by  using  the 
mean  value  of  the  maximum  and  minimum  boundary  of  the  flame 
front.  After  the  flame  made  contact  with  the  wall,  the  geometry 
of  the  flame  deviated  further  from  a  circle.  However,  for  the 
purpose  of  determining  the  center  and  the  radius  of  the  flame, 
the  circle  was  approximated  by  the  boundaries  of  the  flame, 
although  no  attempt  was  made  to  consider  analytically  the 
development  of  the  flame  after  it  had  made  contact  with  the 
chamber  wall. 

Based  on  the  position  of  the  instantaneous  centers  of  the 
flame  thus  determined,  the  flame  trajectory  could  be  obtained 
on  a  reference  rectangular  coordinate  system,  whose  origin 
coincided  with  the  geometrical  center  of  the  chamber,  by 
expressing  each  flame  center  as  a  position  vector  originating 
from  the  center  of  the  chamber.  It  was  possible  to  superimpose 
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4.  9.  FLAME  PATTERN  USED  FOR  THE  DETERMINATION  OF 
FLAME  TRAJECTORY  BY  EXPERIMENTAL  METHOD  AT 
900  RPM.  (3500  FRAMES  PER  SEC.)  (CONTINUOUS 
SINGLE  FILM  STRIP  DISPLAYED  IN  EIGHT  SECTIONS, 

A  THROUGH  H. ) 


FIG.  4.  10  FLAME  PATTERN  WHICH  WAS  USED  TO  FIND  FLAME 
TRAJECTORY  BY  EXPERIMENTAL  METHOD  AT  1,800 
RPM.  (3500  FRAMES  PER  SEC.)  (CONTINUOUS  SINGLE 
FILM  STRIP  DISPLAYED  IN  FOUR  SECTIONS,  A  THROUGH 
D.) 
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the  reference  coordinate  system  on  each  frame  of  the  film 
identically  by  making  use  of  the  perforation  at  the  four  corners 
of  each  frame.  Judging  from  the  extent  of  distortion  in  the 
circular  shape  of  the  flame  shown  in  Fig. 4.  9  and  Fig. 4. 10,  it 
is  believed  that  the  flame  made  contact  with  the  chamber  wall 
at  the  time  represented  by  the  frames  of  Row  £  Column  7  and 
Row  A  Column  4  respectively.  The  triangular  marks  in  Fig.  4.  11 
represent  the  position  of  flame  centers  observed  on  every  other 
frame  of  film  or  at  intervals  of  1.  1  milliseconds.  The  triangular 
marks  in  Fig.  4.  12  represent  the  position  of  flame  centers 
observed  on  every  frame  of  film,  or  at  0.  64  millisecond  intervals 
at  the  camera  speed  used  here.  The  measured  data  of  the  flame 
radius  for  both  speeds  of  combustion  chamber  rotation  are  given 
in  Fig.  4.  13.  The  data  were  obtained  for  the  period  before  the 
flame  made  contact  with  the  chamber  wall. 

Photographs  were  also  taken  for  the  nonhomo geneous  case 
at  a  chamber  speed  of  900  rpm  with  both  the  mixture  charged  and 
ignition  effected  at  half  radius  from  the  center  of  the  chamber. 
These  photographs  are  shown  in  Fig.  4.  14.  As  can  be  seen  in  this 
figure,  a  circular  flame  pattern  did  not  develop. 

4.  6  Results  of  Theoretical  Calculations 

As  noted  in  Equation  2.  9  and  2.  10  given  in  Section  2.  3, 
the  radius  of  the  flame  and  the  value  of  —  as  a  function 

of  time  must  be  known  in  order  to  calculate  the  flame  trajectory. 

Since  both  the  radius  of  the  flame  and  the  combustion 
pressure  are  functions  of  time,  it  should  be  possible  to  determine 
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the  relationship  between  the  radius  of  the  flame  and  the  com¬ 
bustion  pressure.  However,  the  functional  form  of  neither  the 
radius  of  the  flame  nor  the  combustion  pressure  as  a  function 
of  time  is  known  as  yet.  Therefore,  either  the  relationship 
between  the  radius  of  the  flame  and  combustion  time  or  the 
relationship  between  pressure  and  time  has  to  be  determined 
experimentally.  In  the  present  investigation,  it  was  decided  to 
obtain  the  relationship  between  the  radius  of  the  flame  and  the 
combustion  pressure  by  using  Equation  2.  11  derived  by  Lewis 
and  Von  Elbe  (13)  and  the  experimentally  determined  com¬ 
bustion  pressure  as  a  function  of  time.  The  flame  radii 
calculated  in  this  manner  together  with  the  directly  measured 
flame  radii  are  shown  in  Fig.  4.  13.  Since  the  flame  front 
compresses  the  surrounding  acetylene -air  mixture  adiabatic  ally , 
the  temperature,  pressure,  and  density  of  the  mixture  can  be 
found.  Using  these  data,  the  viscosity  coefficient  can  be  found 
from  the  air  tables  under  the  assumption  that  the  viscosity,  of 
the  mixture  is  the  same  as  that  of  air.  Values  of  obtained 

/A 

in  this  manner,  as  a  function  of  time,  are  shown  in  Fig.  4.  15. 
Now  from  these  values  of  the  Reynolds  number  Re  as  a 

function  of  time  during  combustion  was  obtained  from  the 
relationship  shown  in  Section  2.3,  Then  the  drag  coefficient 
was  obtained  from  Re  using  the  linear -logarithmic  approximation 
given  in  Table  4.  1  and  Fig. 4.  16.  With  these  values  for  flame 
radius  R  and  drag  coefficient  (the  calculations  for  both 
have  made  use  of  experimentally  determined  data),  the 
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FIG.  4.  13.  COMPARISON  OF  THE  FLAME  RADIUS 
DETERMINED  BY  THE  THEORY  AND 
EXPERIMENT. 


FIG  4.  14.  FLAME  PATTERN  IN  NONHOMOGENEOUS  MIXTURE 

COMBUSTION  AT  900  RPM.  3500  FRAMES  PER  SECOND. 
0.3  IN.  HG  (ACETYLENE)  WAS  CHARGED  TO  THE 
CHAMBER  MODEL  (PREMIXING  TANK:  21  CU.  IN.) 
ANILINE  WAS  ADDED  TO  THE  PREMIXING  TANK. 
(CONTINUOUS  SINGLE  FILM  STREP  DISPLAYED  IN 
SIX  SECTIONS,  A  THROUGH  F. ) 


FIG.  4. 
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differential  equations  for  the  trajectory  of  the  flame  (Equations 
2.  9  and  2.  10)  were  solved.  The  calculations  were  made  as 
described  in  Appendix  II,  using  an  I.  B.M.  7074  computer  for 
the  theoretical  trajectories  of  the  flame  for  the  three  cases  of 
the  homogeneous  mixture  ignited  at  half  and  three  quarter  radius 
from  the  center  and  at  the  periphery  under  the  combustion  chamber 
rotative  speeds  of  900  and  1800  rpm.  The  results  are  shown  in 
Fig.  4.  17  (fixed  coordinates)  and  Fig.  4.  18  (rotational  coordinates). 
The  results  of  theoretically  and  experimentally  obtained  flame 
trajectories  are  compared  in  Fig.  4.  11  and  Fig.  4.  12. 

For  the  purpose  of  studying  the  effects  of  combustion 
chamber  rotative  speed  and  the  location  of  the  point  of  ignition 
on  the  magnitudes  of  the  position  vector  and  the  tangential  velocity 
of  the  flame  centers,  calculations  were  also  made  as  a  function  of 
time  for  the  above  mentioned  parameters  for  flame  trajectories 
resulting  from  ignition  effected  at  half  and  three-quarter  radius 
from  the  center  and  at  the  periphery  under  the  two  different 
chamber  rotative  speeds  of  900  and  1800  rpm.  The  results  are 
shown  in  Fig.  4.  17  through  Fig.  4.  22,  for  both  fixed  and  rotational 
coordinate  systems.  All  of  these  calculations  were  made  under 
the  assumption  that  the  variation  of  combustion  pressure  with 
time  was  identical  in  all  three  cases  for  the  different  ignition 
points  at  each  of  the  combustion  chamber  rotative  speeds. 

The  results  given  in  this  chapter  will  be  discussed  in 


Chapter  V. 
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Ignition  Position 
Combustion  End 


FIG.  4.  17.  THEORETICAL  FLAME  TRAJECTORY  FOR 
FIXED  COORDINATES. 
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FIG.  4. 18.  THEORETICAL  FLAME  TRAJECTORY  FOR 
ROTATIONAL  COORDINATES. 


Magnitude  of  Position  Vector  of  the  Flame  (in.) 
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FIG.  4.19.  RELATIONSHIP  BETWEEN  THEORETICAL,  AND 
EXPERIMENTAL  MAGNITUDE  OF  POSITION 
VECTOR  OF  THE  FLAME  CENTER  AND  COM¬ 
BUSTION  TIME  FOR  FIXED  COORDINATES. 
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Time  (10"^  Sec) 


FIG.  4.20.  RELATIONSHIP  BETWEEN  MAGNITUDE  OF 
POSITION  VECTOR  OF  THE  FLAME  CENTER 
AND  COMBUSTION  TIME  FOR  ROTATIONAL 
COORDINATES. 


FIG.  4.21.  RELATIONSHIP  BETWEEN  MAGNITUDE  OF 
TANGENTIAL  VELOCITY  AND  COMBUSTION 
TIME  FOR  FIXED  COORDINATES. 
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FIG.  4.22. 


RELATIONSHIP  BETWEEN  MAGNITUDE  OF 
TANGENTIAL  VELOCITY  AND  COMBUSTION 
TIME  FOR  ROTATIONAL  COORDINATES. 
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V.  DISCUSSION  OF  EXPERIMENTAL  AND 

- THEORETICAL  RESULTS - 

5.  1  Discussion  of  Experimental  Results 

When  a  flame  occurs  in  a  swirling  field,  the  flame  moves 
inward  to  the  center  of  the  chamber  due  to  an  unbalance  of  the 
pressure  and  drag  forces.  These  forces  are  considered  here  to 
be  the  only  significant  forces  in  the  system. 

As  the  flame  moves  inward  during  the  combustion,  the  air- 
fuel  mixture  will  displace  outward  due  to  the  pressure  force,  and 
this  motion  produces  turbulence  in  the  chamber.  This  turbulence 
increases  the  mixing  effect  between  air  and  fuel  even  in  the  case 
of  a  mixture  which  was  originally  homogeneous  before  ignition 
has  occurred. 

Flame  motion  changes  with  the  change  of  rotative  chamber 
speed  since  the  latter  results  in  a  change  of  pressure  and  drag 
forces.  Furthermore,  as  the  combustion  time,  flame  motion, 
as  related  to  the  combustion  time,  is  dependent  upon  the 
chamber  rotative  speed.  Figure  4.  1  clearly  indicates  that  the 
pressure  pattern  is  considerably  influenced  by  rotative  chamber 
speed.  These  data  imply  that  the  combustion  time  is  greatly 
influenced  by  the  combustion  chamber  rotative  speed,  in  view  of 
the  definition  of  combustion  time  given  earlier.  Throughout  these 
test 8,  the  maximum  combustion  pressure  was  found  to  be  nearly 
constant  when  the  same  amount  of  mixture  was  charged.  That 
is,  the  loss  of  energy  is  nearly  constant  during  the  combustion  at 
different  speeds. 
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As  shown  in  Fig. 4.  5,  when  ignition  takes  place  at  the 
periphery,  the  combustion  time  is  influenced  by  the  chamber 
rotative  speeds.  The  pressure  force  acting  on  the  flame  in  the 
case  of  high  speeds  Is  larger  than  that  of  the  case  of  stationary 
and  low  chamber  rotative  speeds;  thus,  showing  that  the  pressure 
force  has  a  definite  effect  on  combustion  time.  At  comparatively 
high  speeds,  the  change  in  combustion  time  becomes  small  since 
the  flame  front  moves  very  rapidly  due  to  the  greater  pressure 
force,  without  experiencing  perfect  combustion.  The  gas  which 
could  not  burn  perfectly  is  converted  to  smoke  behind  the  flame 
front  and  therefore  the  change  in  combustion  time  with  respect  to 
combustion -chamber  rotative  speed  becomes  small  at  high  speeds. 
The  difference  between  the  combustion  time  for  peripheral,  half 
radius,  and  center  ignition  as  indicated  in  Fig.  4.  5,  is  mainly  due 
to  the  differences  in  the  area  of  the  burning  surface  and  in  the  wall 
effect.  This  is  true  both  in  the  case  of  homogeneous  and  nonhomo- 
geneous  mixture  combustion.  In  the  case  of  peripheral  ignition,  the 
flame  configuration  is  not  cylindrical  and  the  burning  surface 
area  is  smaller  than  in  the  cases  of  half  radius  from  the  center 
or  center  ignitions.  The  time  variation  of  the  pressure,  which 
was  recorded  with  the  Polaroid  camera,  is  not  much  different  until 
the  wall  effect  becomes  important,  as  shown  in  Fig.  5.  1.  The  graph 
was  obtained  for  methane -air  mixture  charged  into  the  stationary 
chamber  in  the  case  of  Part  (c).  In  the  case  of  center  ignition 
(refer  to  the  description  in  3.  1.  1),  the  chamber  rotative  speed 
does  not  influence  the  combustion  time  because  the  pressure 
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FIG.  5. 1.  COMPARISON  OF  COMBUSTION  PRESSURE  VS. 

COMBUSTION  TIME  WHEN  IT  WAS  IGNITED  AT 
HALF  RADIUS  FROM  CENTER  AND  CENTER. 
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is  zero  at  the  geometric  center.  Therefore  a  trajectory  does  not 
exist  in  this  case,  and  the  flame  will  propagate  radially.  This  is 
equivalent  to  the  case  of  the  fundamental  Equation  2.  9  with  all 
the  boundary  conditions  equal  to  zero,  which  makes  the  case 
trivial. 

As  shown  in  Fig.  4.  6,  as  the  chamber  rotative  speed  is 
increased,  the  combustion  time  decreases  for  each  of  the  three 
ignition  positions  with  nonhomo geneous  mixture  combustion  and 
with  the  mixture  charged  at  the  periphery.  The  reason  is  that 
during  combustion,  the  nonhomogcneous,  over  rich  mixture  at 
high  chamber  speeds  becomes  a  comparatively  lean  mixture  in 
the  central  area  of  the  chamber  due  to  forced  diffusion.  It  is 
believed  that  such  a  lean  mixture  resulting  at  high  speeds  is 
more  suitable  for  combustion  than  the  richer  mixture  present  at 
low  speeds.  The  speed  of  the  chamber  has  a  great  effect  on  the 
combustion  time  in  the  case  of  center  ignition;  for,  as  the  rich 
mixture  was  charged  at  the  periphery,  diffusion  occurred  due 
to  the  pressure  gradient  which  has  a  direct  relationship  to  the 
chamber  speed.  The  mixture  at  the  geometrical  center  becomes 
lean  compared  with  the  mixture  around  the  charging  port. 

This  condition  is  believed  to  be  more  suitable  for  combustion; 
that  is,  the  flame  which  originates  at  the  geometrical  center 
will  propagage  itself  radially  with  a  high  speed  at  the  prevailing 
lean  mixture  ratio. 

When  the  fuel  was  charged  at  half  radius  from  the  center, 
the  effect  of  combustion  chamber  rotative  speed  on  combustion 
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time  was  exactly  opposite  from  that  for  the  results  shown  in 
Fig.  4.  6,  obtained  when  charging  was  effected  at  the  periphery. 

As  shown  in  Fig.  4.  7,  combustion  time  increases  with  an  increase 
in  chamber  speed.  The  reason  for  this  is  that  part  of  the  rich 
mixture  diffuses  toward  the  geometrical  center,  due  to  the 
pressure  gradient,  while  another  portion  of  the  rich  mixture 
will  diffuse  toward  the  wall  due  to  the  temperature  gradient 
caused  by  the  heat  of  friction  between  the  journal  bearing  and 
hollow  cylinder  produced  by  the  chamber  rotation.  But  the 
main  portion  of  diffusion  is  believed  to  be  the  forced  diffusion 
due  to  the  pressure  gradient  rather  than  the  thermal  diffusion. 
When  the  point  of  ignition  was  located  at  the  periphery  of  the 
chamber,  the  flame  around  the  periphery  propagated  more 
rapidly  than  that  around  the  charging  port;  and  as  the  flame 
front  arrived  near  the  center  where  acetylene  was  concentrated 
due  to  force  diffusion  from  the  charging  port,  flame  speed 
decreased  further  from  that  around  the  periphery.  Finally, 
the  flame  passed  through  the  center  and  propagated  toward  the 
lean  mixture  zone  formed  near  the  wall  opposite  to  the  ignition 
point.  Thus  the  most  unsuitable  region  for  combustion  was 
considered  to  be  around  the  center  in  the  case  of  charging  at 
half  radius  from  the  center.  In  other  words,  the  concentration 
of  the  acetylene-air  mixture  around  the  center  increased  with 
the  increase  in  chamber  speed  and  caused  the  flame  speed  to 
decrease  around  the  center. 
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Therefore  combustion  time  increased  with  the  increase  of 
rotative  chamber  speed.  For  a  given  value  of  rotative  chamber 
speed,  the  combustion  time  was  found  to  be  the  shortest  for  the 
center  ignition  and  the  largest  for  the  peripheral  ignition.  This 
fact  may  be  due  to  the  larger  burning  surface  area  that  occurs 
with  center  ignition  and  the  smaller  burning  surface  area  that 
occurs  with  peripheral  ignition. 

As  shown  in  Fig. 4.  8,  in  the  case  of  homogeneous  mixture 
combustion  of  a  rich  mixture,  the  chamber  rotative  speed  has 
more  effect  on  combustion  time  than  that  in  the  case  of  a  non- 
homogeneous  mixture  combustion.  This  condition  is  due  to  the 
fact  that,  in  the  case  of  homogeneous  mixture  combustion,  the 
higher  the  combustion  chamber  rotative  speed  the  higher  the 
pressure  gradient  becomes,  and  in  turn,  the  higher  the  speed 
of  flame  propagation  becomes. 

In  the  case  of  nonhomogene ous  mixture  combustion  the 
mixture  ratio  varies  throughout  the  chamber  due  to  the  pressure 
gradient  and  the  temperature  gradient.  The  distribution  of  the 
mixture  ratio  in  the  nonhomo geneous  case  is  better  for  com¬ 
bustion  than  that  of  the  homogeneous  mixture.  Based  on  the  above 
given  argument,  it  can  be  said  that  the  relationship  between  com¬ 
bustion  time  and  chamber  rotative  speed  is  less  pronounced  in  the 
case  of  nonhomo  gene  ous  mixture  combustion  than  in  the  case  of 
homogeneous  mixture  combustion. 
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As  a  counterpart  of  the  experiments  carried  out  with  an 
overrich  mixture  Part  (e)  and  (g),  a  similar  experiment  was 
carried  out  with  a  lean  acetylene -air  mixture  with  a  ratio  of 
1:6.4  (acetylene  of  10  inches  of  mercury  above  atmospheric 
pressure  and  air  of  28  inches  of  mercury  above  atmospheric 
pressure)  charged  into  the  combustion  chamber  (0.3  x  10  lb). 
Combustion  time  was  found  to  be  shorter  for  the  homogeneous 
mixture  combustion  than  for  the  nonhomogeneous  mixture  com¬ 
bustion.  In  the  case  of  nonhomogeneous  mixture  combustion, 
combustion  time  increased  with  increasing  chamber  speed; 
whereas,  in  the  homogeneous  mixture  combustion,  combustion 
time  decreased  with  an  increase  in  the  rotative  chamber  speed. 
The  reason  can  be  explained  as  follows.  A  mixture  ratio  of 
1:6£  is  close  to  a  stoichiometric  ratio;  and,  in  the  case  of 
homogeneous  mixture  combustion,  combustion  time  decreased 
with  an  increase  in  chamber  rotative  speed  due  to  the  pressure 
force.  However,  in  the  case  of  nonhomogeneous  mixture 
combustion,  the  mixture  ratio  became  leaner  toward  the  center 
due  to  forced  diffusion;  and  near  the  center,  the  mixture  ratio 
decreased  to  the  extent  that  the  speed  of  flame  propagation  is 
reduced  greatly  (14). 

The  above  explanation  leads  to  the  conclusion  that  non¬ 
homogeneous  mixture  combustion,  such  as  that  in  a  high  speed 
diesel  engine,  can  provide  larger  values  of  time  rate  of  energy 
emission  than  the  homogeneous  mixture  combustion,  provided 
an  over  rich  mixture  is  used  at  high  swirl  rates. 
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In  the  experiment  with  nonhomogeneous  mixture  com¬ 
bustion,  a  pressure  spike  was  observed  during  the  combustion 
process.  In  order  to  determine  the  origin  of  this  spike,  a  one 
dimensional  cylindrical  chamber  (with  a  length  much  greater 
than  the  diameter),  shown  in  Fig.  5.  2,  was  constructed.  First 
a  pressure  pickup  with  a  connecting  tube  was  attached,  and  then 
a  pressure  pckup  without  a  connecting  tube  was  attached  to  this 
combustion  chamber  in  an  effort  to  determine  the  cause  of  the 
pressure  spike.  The  experiment  was  carried  out  for  the  non¬ 
homogeneous  case  using  both  rich  and  lean  mixtures.  The  rich 
and  lean  acetylene-to-air  mixture  ratio  (by  weight)  was  1:4,85 
and  1:6.4  respectively.  Fuel  and  air  were  mixed  in  the  pre¬ 
mixing  tank  and  charged  into  the  combustion  chamber  until  the 

pressure  in  the  premixing  tank  decreased  to  2  inches  of 

-4 

mercury  equivalent  to  an  outflow  of  0.  605  x  10  lb.  The 
relationship  between  combustion  pressure  and  time  was  recorded 
using  the  Polaroid  camera. 

The  results  are  shown  in  Fig.  5.  3  (a)  and  (b) ,  where  it 
is  noted  that  a  pressure  spike  was  observed  only  during  the 
tests  with  the  pressure  pickup  which  had  a  connecting  tube. 

The  results  of  these  tests  indicate  the  following  to  be  the 
cause  of  the  pressure  spike.  The  small  tube  which  is  connected 
to  the  pressure  pickup  will  be  filled  with  combustible  gas  when 
an  over  rich  mixture  is  charged  into  the  chamber  and  it  will 
burn  out  as  the  flame  front  passes  the  geometrical  center  of 
the  chamber.  The  burning  gas  will  cause  the  pressure  in  the 


FIG.  5.  2.  ONE  DIMENSIONAL  COMBUSTION  CHAMBER  MODEL  FOR  THE  PRESSURE 

SPIKE. 
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FIG.  5.  3.  PHOTOGRAPHIC  RECORD  OF  PRESSURE  SPIKE  PHENOMENA. 
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tube  to  rise  very  rapidly  because  the  volume  of  the  tube  is  small. 
The  burning  gas,  with  rising  pressure,  will  produce  a  jet  into 
the  main  combustion  chamber  due  to  the  pressure  difference 
between  the  tube  and  the  main  combustion  chamber  and  cause  the 
pressure  to  decrease  very  rapidly  in  the  tube.  Thus  it  can  be 
concluded  that  the  occurence  of  the  pressure  spike  is  due  to  a 
phenomenon  that  has  no  effect  on  the  combustion  time  in  cases 
as  (d)  of  Fig.  5.  3.  The  pressure  spike  did  not  appear  with  either 
of  the  pressure  pickup  installation  in  lean  mixture  combustion 
because  the  mixture  in  the  tube  did  not  burn. 

5.  2  Discussion  of  Photographic  Results 

As  shown  in  Fig.  4.  2  and  Fig.  4.  3  for  the  homogeneous 
case,  a  strong  swirl  was  observed  as  the  flame  disappeared. 

This  swirl  is  considered  to  be  caused  by  the  inward  move¬ 
ment  of  the  flame  during  the  early  stages  of  the  combustion. 

When  the  acetylene -air  mixture  as  specified  in  Part  (h)  of 
Table  3.  1  was  charged  into  the  combustion  chamber  at  the  half 
radius  from  the  center  and  ignited  at  half  radius  at  900  rpm,  the 
flame  was  observed  to  spiral  around  the  geometrical  center  in 
the  later  stage  of  combustion.  The  experimental  methods  are 
not  idealized  enough  to  permit  comparison  with  the  theoretical 
results  which  were  obtained  under  the  assumption  that  the  flame 
is  circular. 

In  finding  the  flame  trajectory  experimentally,  it  is 
advantageous  to  study  the  nonhomogeneous  mixture  combustion  in 
which  the  flame  cam  propagate  symmetrically  about  the  flame 


88 


center.  In  the  nonhomogeneous  case,  the  rich  mixture  slug 
charged  into  the  chamber  will  remain  circular  in  the  ideal  case, 
and  one  can  expect  to  have  a  flame  which  expands  adiabatic  ally 
and  symmetrically  about  the  flame  center  in  the  mixture.  Then, 
as  the  flame  expands,  the  boundary  of  the  acetylene -air  mixture 
will  be  displaced  outward  by  the  flame  front.  The  flame  center 
will  deviate  from  the  initial  center  due  to  the  pressure  force. 

At  the  last  stage  of  combustion,  the  flame  front  will  coincide 
with  the  boundary  of  the  expanding  mixture.  In  the  diesel  case, 
the  form  of  the  flame  is  circular;  and  this  flame  will  move  for 
a  while  until  it  diminishes.  Consequently,  one  can  observe 
the  flame  trajectory  by  photographing  the  flame  with  a  high 
speed  camera.  In  this  case,  if  the  maximum  radius  of  the 
flame  is  small  compared  to  the  radius  of  the  combustion  chamber, 
the  drag  force  will  be  proportional  to  the  square  of  the  relative 
velocity  of  the  flame.  Thus  the  original  differential  equation 
will  apply,  even  for  the  later  portion  of  the  combustion.  This 
condition  may  be  true  ideally;  however,  in  the  experimental 
work  it  was  found  that,  for  the  nonhomogeneous  case,  the  flame 
was  not  quite  circular,  as  shown  in  Fig.  4.  14.  This  caused  the 
measured  results  to  differ  from  the  theoretical  predictions. 
Therefore,  the  measurements  for  the  homogeneous  case, 
ignited  at  the  point  half  radius  from  the  center,  shown  in  Fig.  4.  9 
and  Fig,  4.  10,  were  used  for  comparison  with  the  theoretical 
predictions.  It  can  be  seen  that  in  the  early  stage  of  the  com¬ 
bustion  the  flame  follows  the  original  assumption  that  the  flame 
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is  circular  until  the  wall  effect  becomes  appreciable, 

5.  3  Discussion  of  Theoretical  Results 

As  shown  in  the  Fig.  4.  16  and  Fig.  4.  17,  the  trajectory 
of  the  flame  in  the  case  of  high  chamber  rotative  speeds  is  almost 
the  same  as  that  of  the  case  of  lower  speeds.  If  the  drag  force 
is  neglected  and  the  pressure  force  is  considered  as  the  only 
force  acting  on  the  system,  the  trajectory  of  the  flame  has  no 
relation  to  the  chamber  speed  as  shown  by  Equation  2.  15. 
Conversely,  if  we  neglect  the  pressure  force  and  consider  the 
drag  force  as  the  only  force  acting  on  the  flame,  the  flame 
trajectory  (solution  of  Equation  2.  16)  is  found  to  be  almost 
independent  of  the  combustion  chamber  rotative  speed.  There¬ 
fore,  we  can  postulate  that  the  actual  flame  trajectory  is  almost 
independent  of  the  chamber  rotative  speed.  However,  the  rate 
of  change  of  the  radius  of  curvature  of  the  flame  trajectory  in 
the  case  of  1800  rpm  is  almost  twice  that  of  the  case  of  900  rpm. 
(Refer  to  Fig.  4,  19  and  Fig.  4. 20) 

If  we  consider  that  both  the  pressure  and  drag  forces  are 
acting  on  the  flame,  the  trajectory  of  the  flame  will  deviate  from 
an  elliptical  curve  toward  the  inside  due  to  the  drag  force  as 
shown  in  Fig.  5.4  (solution  of  Equation  2.  9) .  The  deviation  of 
the  trajectory  from  an  elliptical  path  is  very  slight  at  the 
beginning  of  combustion  since  the  drag  force  acting  on  the  flame 
is  very  small  because  the  tangential  velocity  is  almost  the  same 
as  the  rotative  chamber  speed  at  the  beginning  of  combustion, 
as  shown  in  Fig.  5.  4  and  Fig.  5.  5. 
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FIG.  5.  4.  COMPARISON  OF  FLAME  TRAJECTORY  IN  CASE 
OF  ONLY  PRESSURE  FORCE  ACTING  OR  ONLY 
DRAG  FORCE  ACTING  AS  THE  FICTITIOUS  CASE 
AND  COMBINED  FORCES  ACTING  ON  THE  FLAME 
FOR  FIXED  COORDINATE  SYSTEM. 
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COMBUSTION  TIME  IN  CASE  OF  ONLY  PRESSURE 
FORCE  ACTING  OR  ONLY  DRAG  FORCE  ACTING 
AND  COMBINED  FORCES  ACTING  ON  THE  FLAME. 
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The  radius  of  curvature  of  the  flame  trajectory  becomes 
smaller  as  the  mixture -to -flame  density  ratio  increases,  as  shown 
in  Fig.  4.  17  and  Fig.  4.  18.  This  change  is  due  to  the  unbalanced 
forces.  This  conclusion  was  reached  under  the  assumption  that 
the  combustion  time  and  pressure  are  the  same  as  that  in  the  case 
of  ignition  at  the  point  half  a  radius  from  the  center.  As  shown 
in  Fig.  4.  20,  the  magnitude  of  the  position  vector  of  the  flame  is 

S<L 


smaller  in  the  case  where 


=  4.0  than  that  of  the  case 


where 


3.0. 


The  tangential  velocity  of  the  flame  for  the  combustion 
chamber  rotative  speed  of  1800  rpm  is  almost  twice  that  of  the 
case  when  the  chamber  rotative  speed  was  900  rpm,  as  shown 
in  Fig.  4.  21  and  Fig. 4.  22.  This  means  that  the  drag  force, 
which  affects  the  mixing  process,  will  increase  by  a  factor  of 
about  four  because  the  drag  force  is  proportional  to  the  square 
of  relative  velocity.  Furthermore,  the  results  shown  in  Fig.  4.  21 
and  Fig.  4.  22  imply  that  stronger  mixing  of  the  fuel-air  mixture 
is  produced  if  ignition  is  effected  at  a  larger  radial  distance  from 
the  center.  This  conditions  is  so  because  fuel-air  mixing  depends 
upon  the  tangential  velocity  of  the  flame  which  is  affected  by  the 
position  of  the  point  of  ignition. 

The  tangential  velocity  increases  as  the  mixture  to  flame 
density  increases  as  shown  in  Fig.  4.  22. 

If  we  neglect  the  drag  force  and  consider  the  hypothetical 
case  where  the  pressure  force  is  the  only  force  acting  on  the 
flame,  then  the  absolute  value  of  the  position  vector  of  the  flame 
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center  is  periodic  and  will  have  a  minimum  value  at  (0, 
0.625/1/3-^  as  shown  in  Fig.  5.4.  In  the  practical  case --that 
is,  when  both  pressure  and  drag  forces  are  acting--the  magnitude 
of  the  position  vector  has  a  minimum  value  that  is  smaller  than 
that  of  the  above-mentioned  hypothetical  case  because  of  the  drag 
force.  The  drag  force  corresponds  to  a  damping  force  in  a  mass- 
spring  system. 

As  observed  in  Equation  2.  15,  when  the  mixture  -to  -flame 
density  ratio  increases,  the  flame  center  describes  a  slender 


ellipse.  For  example,  when  ignition  occurs  at  (0.625,  0)  with 

— ^ —  =4.0,  the  major  and  minor  axes  of  the  ellipse  are 

Jf 

0.625  and  0.625/2  respectively.  In  the  ultimate  case,  if  the 
Jcls 

density  ratio  (-5—)  becomes  infinite,  the  flame  motion  becomes 
a  simple  harmonic  motion.  Therefore,  one  can  predict  that  a 
flame  occurring  in  a  mixture  which  has  larger  value  of  will 

move  inside  of  the  trajectory  of  a  flame  occurring  in  a  mixture 
having  a  smaller  value  of 

As  shown  in  Fig.  4.  22,  the  magnitude  of  the  tangential 


=  4.0,  is  larger  than  that  of  the  case  of 


velocity  of  the  flame,  in  a  rotational  coordinate  system,  and  in 
Scu 

the  case  of  — 3 —  =4.0,  is  larger  than  that  of  the  case  of 

A.  fs 

— ; 3 —  =  3.0.  Therefore,  in  the  practical  case,  it  can  be 

predicted  that  combustion,  where  — s =4.0,  is  better 

JL  * 

than  that  where  — s —  =3.0  because  the  greater  drag  force 


predicted  that  combustion,  where 

i.1 _ J-L.i,  _  AT _ 


=  4.0,  is  better 


acts  in  the  former  case,  and  the  greater  drag  force  causes 


better  mixing  and  therefore  better  combustion. 

According  to  Fig.  5.4,  the  flame  trajectory  resulting  from 
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the  pressure  force  alone  almost  coincides  with  the  flame  trajectory 
obtained  when  the  combined  forces  are  considered;  whereas,  the 
trajectory  resulting  from  the  drag  force  alone  deviates  considerably 
from  that  of  the  combined  forces.  This  leads  to  the  conclusion 
that  the  flame  trajectory  is  determined  predominantly  by  the 
pressure  force. 

As  a  hypothetical  case,  if  one  considers  that  the  pressure 
force  is  the  only  force  acting  on  the  flame,  the  relationship 
between  the  tangential  velocity  of  the  flame  and  time  is  a  periodic 
function  of  time,  and  the  flame  motion  follows  Kepler’s  second 
principle,  as  shown  in  Fig.  5..  5. 

In  the  case  where  the  drag  force  is  considered  as  the  only 
force  acting  on  the  flame,  the  tangential  velocity  is  0.625 
which  is  comparable  to  the  combustion  chamber  rotative  speed  at 
the  beginning  of  combustion;  that  is,  at  the  beginning  the  drag 
force  is  almost  negligible,  as  shown  in  Fig.  5.  5.  As  the  drag 
force  increases  withtime,  the  flame  deviates  from  the  tangential 
line  at  (0.  625,  0) ,  and  the  rate  of  deviation  from  the  tangential 
line  in  Fig.  5.4  will  cause  the  actual  flame  trajectory  to  describe 
a  spiral  curve  inside  of  the  elliptical  trajectory  which  was  obtained 
from  the  hypothetical  case. 

As  shown  in  Fig.  5.5,  the  magnitude  of  the  tangential 
velocity  of  the  flame  will  increase  to  a  maximum  which  is  caused 
by  the  pressure  force.  As  the  drag  force  increases,  this  increase 
will  decelerate  the  magnitude  of  the  tangential  velocity  of  the 


flame . 
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5.  4  The  Comparison  of  the  Experimental  and  Theoretical  Results 
for  the  Flame  Trajectory. 

Figure  4.  17  and  Fig.  4.  18  show  the  flame  trajectory  which 
was  obtained  by  theoretical  calculations  for  fixed  and  rotational 
coordinates.  Figure  4.11  and  Fig.  4.  12  show  the  comparison  of 
the  theoretical  and  experimental  trajectories  of  the  flame.  At 
the  beginning  of  combustion,  both  curves  coincide;  but,  in  the 
later  portion  of  the  combustion,  the  curves  deviate  considerably 
due  to  the  wall  effect  and  the  drag  force.  In  the  calcaulation  of 
the  drag  force,  it  was  assumed  that  the  boundary  line  of  the  flame 
is  distinct;  but,  in  the  practical  case  the  boundary  line  is  not 
distinct;  and,  when  the  flame  grows  to  a  large  size  compared  with 
the  size  of  the  chamber,  the  drag  force  is  no  longer  proportional 
to  the  square  of  the  relative  velocity.  The  expanding  flame  does 
not  follow  the  adiabatic  process,  see  Fig.  4.  14.  As  shown  in 
Fig.  5.  5,  the  experimental  magnitude  of  the  tangential  flame 
velocity  is  smaller  than  that  of  the  theoretical  calculation;  also, 
as  shown  in  Fig.  4.  20,  the  magnitude  of  the  measured  position 
vector  of  the  flame  is  smaller  than  that  of  the  theoretical 
calculation  because  the  drag  force  acting  in  the  practical  case 
is  larger  than  that  assumed  in  the  theoretical  case,  even  in 
the  early  stages  of  the  combustion. 
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coordinates.  Figure  4.11  and  Fig.  4.  12  show  the  comparison  of 
the  theoretical  and  experimental  trajectories  of  the  flame.  At 
the  beginning  of  combustion,  both  curves  coincide;  but,  in  the 
later  portion  of  the  combustion,  the  curves  deviate  considerably 
due  to  the  wall  effect  and  the  drag  force.  In  the  calcaulation  of 
the  drag  force,  it  was  assumed  that  the  boundary  line  of  the  flame 
is  distinct;  but,  in  the  practical  case  the  boundary  line  is  not 
distinct;  and,  when  the  flame  grows  to  a  large  size  compared  with 
the  size  of  the  chamber,  the  drag  force  is  no  longer  proportional 
to  the  square  of  the  relative  velocity.  The  expanding  flame  does 
not  follow  the  adiabatic  process,  see  Fig. 4.  14.  As  shown  in 
Fig.  5.  5,  the  experimental  magnitude  of  the  tangential  flame 
velocity  is  smaller  than  that  of  the  theoretical  calculation;  also, 
as  shown  in  Fig.  4.  20,  the  magnitude  of  the  measured  position 
vector  of  the  flame  is  smaller  than  that  of  the  theoretical 
calculation  because  the  drag  force  acting  in  the  practical  case 
is  larger  than  that  assumed  in  the  theoretical  case,  even  in 
the  early  stages  of  the  combustion. 
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VI,  SUMMARY  AND  CONCLUSION 

6.  1  Statement  of  the  Problem 

The  purpose  of  this  investigation  was  to  investigate  the 
trajectory  of  the  flame  in  a  swirling  field  by  theoretical  and 
experimental  methods.  In  addition,  the  effect  of  the  rotative 
combustion  chamber  speed  on  the  combustion  time  in  homo¬ 
geneous  and  nonhomogeneous  mixture  combustions  was  investi¬ 
gated. 

The  results  of  the  theoretical  work  were  compared  with 
the  experimental  results  in  finding  the  trajectory  of  the  flame  in 
the  swirling  field  of  the  combustion  chamber. 

6.  2  Procedure  and  Results 

The  first  part  of  this  work  was  the  investigation  of  the 
effects  of  the  rotative  combustion  chamber  speed  on  the  com¬ 
bustion  time  in  homogeneous  and  nonhomogeneous  mixtures  using 
both  methane  and  acetylene  gas. 

The  second  part  of  this  work  was  the  investigation  of  the 
trajectory  of  the  flame  in  homogeneous  mixtures  using  acetylene 
gas  as  fuel.  A  high  speed  camera  was  used  to  record  the  visual 
observations. 

Test  results  show  that  the  combustion  chamber  rotative 
speed  has  an  influence  on  the  combustion  time.  It  has  been  found 
that  the  optimum  location  of  the  point  of  ignition  is  at  the  periphery 
in  both  homogeneous  and  nonhomogeneous  mixture  combustion. 

In  the  homogeneous  mixture  combustion,  the  effects  of  combustion 


97 


chamber  rotative  speed  on  the  combustion  time  was  greater  than 
that  of  the  nonhomogeneous  mixture  combustion  for  rich  mixtures, 
but  the  condition  of  combustion  in  the  nonhomogeneous  case  was 
better  than  that  of  the  homogeneous  mixture  combustion.  For 
lean  mixtures,  the  homogeneous  mixture  has  better  combustion 
conditions  than  that  of  the  nonhomogeneous  mixture  combustion. 

The  theoretical  trajectory  of  the  flame  almost  coincides 
with  the  experimentally  determined  trajectory  of  the  flame  at 
the  beginning  portion  of  the  combustion.  If  the  density  ratio  of 
the  air  to  the  flame  increases,  one  can  predict  that  the  combustion 
time  will  decrease  and  thermal  mixing  will  improve. 

Comparing  the  magnitude  of  the  tangential  velocity  of  the 
theoretical  and  experimental  results,  it  was  found  that  the  drag 
force,  in  the  experimental  case,  was  larger  than  that  of  the 
theoretical  calculation. 

The  larger  deviation  between  the  theoretical  and  experi¬ 
mental  flame  trajectories  at  the  later  portion  of  the  combustion 
may  be  caused  by  the  inaccurate  calculation  of  the  drag  force  and 
the  wall  effect. 

6.  3  Conclusions 

The  optimum  ignition  position  is  on  the  periphery  of  the 
chamber  because  the  flame  can  then  introduce  the  strongest  swirl 
in  the  chamber  for  both  the  case  of  homogeneous  and  of  non¬ 
homogeneous  mixture  combustion. 
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It  can  be  expected  that  the  nonhomo  gene  ous  case  will  pro¬ 
vide  better  burning  conditions  than  that  of  the  homogeneous  mixture 
combustion.  The  swirling  effect  in  the  homogeneous  mixture  com¬ 
bustion,  however,  is  better  than  that  of  the  nonhomogeneous  mix¬ 
ture  combustion  for  rich  mixtures. 

At  the  beginning  of  the  combustion,  the  experimentally 
observed  trajectory  of  the  flame  coincides  with  the  theoretical 
trajectory;  thus  the  original  assumptions  can  apply  to  the  practical 
problem  during  the  early  period  of  the  combustion.  But  in  the 
later  portion  of  the  combustion,  the  assumptions  for  the  drag 
force  made  in  this  study  must  be  modified. 

6.  4  Suggestions  for  Further  Work 

It  is  very  important  to  investigate  the  flame  trajectories 
as  they  occur  in  actual  engines  as  a  means  of  improving  com¬ 
bustion.  Some  theoretical  analogy  should  be  found  that  will  more 
nearly  approximate  the  conditions  found  in  actual  engines.  When 
the  flame  grows  to  a  larger  size,  the  drag  force  is  no  longer 
proportional  to  the  square  of  the  relative  velocity.  Thus,  a  number 
of  areas  exist  wherein  a  more  precise  analysis  is  needed  in  order 
to  obtain  a  better  determination  of  the  drag  force. 
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APPENDIX  I 

MATHEMATICAL  DERIVATION 


Numerical  Solution  of  The  Flame  Trajectory  Equations 

From  Chapter  II,  Equation  2*  5,  the  original  equation  of 


the  flame  trajectory  for  fixed  coordinates  is: 


/hi 


Dividing  this  into  its  real  and  imaginary  parts  results  in: 

x t  i  l “■ '  s?1  Vti+iKT+ij h**>f 

..her.  "  3 

3  C 

j (X+frof-t  ( -  3^Ch)2  . 

}=  -  5^\/ci+K+(H"3* 

Applying  the  same  procedure  to  Equation  2,  8  Chapter  II,  for  the 
rotational  coordinate  system  results  in: 


A3 


The  generalized  form  of  Equations  A.  2  and  A.  3  is  as 


follows: 


i=~  f  c  ^3,*. 

}— /OU.-0O 

i=F(ii  t  f) 

$  =<r 


A-*/ 
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From  the  Kutta-Simpson  one-third  rule,  Equations  A,  4 


can  be  expressed  as  a  generalized  set  of  simultaneous  differential 


equations 


cu 


or 

i£=p(p,B  IS) 
P/S  5) 


wJiere. 


Using  the  generalized  Runge-Kutta  method  (15) 

4  p=  (*'p  42  7°  7^*P 
4p=  -£(4pi2.4*p  42.4*P+4,yF) 

*£=  -£(*'*-  +24"jC424//^44'*/X  ) 

4h-r(*y+**7  t**7“*fi 

4'p  =  f(  ?■>) 

4"P  =fW»+i-A>^  p.ti“3'p„,T.+i<:p., 'J.+S^K 
^*p — f  (=(.+i-A*p4,  f.+i'S'B .  P»+i«''p5 

f^+KP-A, R ii*%  P.+i«% *.+i<*8<0i/. 
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^  P  L  Pj  A*-  4^*  are  found  similarly  to  theitfP - &'% 

All  the  4^5  must  be  found  first,  then  all  the  /£*  and  so  on. 
Therefore,  at  t  =  h  +  o 

P=Po+4P  .  +** 

_  -  yj 

p=fi,+*P  W*+4fr 


This  has  been  programmed  for  solution  on  an  I.B.M.  7074 
computer  as  described  in  Appendix  II. 


Calculation  of  the  radius  of  the  flame 

The  notations  used  here  are  defined  as  follows: 


c 

specific  heat  at  constant  pressure 

Btu/lb 

M 

initial  gas  mass  (air -acetylene) 

lb 

m 

mb 

n 

burned  gas  mass 

fraction  of  burned  gas  =  .  o  +  n*l 

radius  of  the  chamber 

lb 

m 

R 

in. 

T 

temperature  *R 

V 

volume  of  the  chamber 

3 

in. 

4 

radius  of  the  flame 

in. 

v4 

volume  of  burned  gas 

3 

in. 

oC  C«/ci, 


(subscripts  i\  and  t>  are  unburned  and  burned  gas 
respectively) 

When  the  process  can  be  assumed  to  be  an  isotropic  process,  as 
in  our  case,  Campell  (9)  derived 
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(fiA  - 


o(~  -ft/"* 


fa  -1 

r* 

CbtA 

Cvi* 


AS 


When  this  equation  is  not  sensitive  to  the  value  of  *0 

p-PV 

which  is  true  in  our  case,  the  simple  equation  'Y)—  "i: - 

Pe  -Pi 

is  correct  to  within  3%  for  practical  purposes. 

The  initial  mass  is  given  by 


M|=  \°.¥- 

1  R  To 


K  I  Cov^rtcl'ht 


and  the  volume  of  the  burned  gas  is  given  by 


K  = 


A-'l 


Substituting  Equation  A.  6  into  Equation  A.  7  we  have 


\/  =  i/  r  i  -  —  —  n 

b  11  To  P  Tte-ftJ 


A.  3 


-—pi-  Zi  f  To  \ 

L  1  Te>  (  P  '  Pe-fe  J 


AA 
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Now  in  an  isentropic  process  we  can  write 


p 

p.  > 


Therefore 


Jl 

Jo 


p_  ^ 

P<>  ' 


Substituting  Equation  A.  10  into  Equation  A.  9,  we  have 
for  the  radius  of  the  flame  as  a  function  of  chamber  pressure 
and  time 


K»— fcO  (-£-)  fy.w 

This  is  Equation  2.  1 1  in  Chapter  II  which  was  used  to  determine 
the  radius  of  the  flame  for  use  in  the  flame  trajectory  equations. 
Equations  2.9,  2.  10,  and  2.16. 
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APPENDIX  II 


COMPUTER  PROGRAM  FOR  EVALUATING 


Since  a  hand  solution  of  the  Equations  2.  9  and  2,  10  in 
Chapter  II  would  have  been  too  laborious,  it  was  decided  to  use 
the  University's  computer,  anl.B.M.  7074. 

The  program  for  the  solution  of  the  equation  by  the 
generalized  Runge-Kutta  method  was  written  in  Fortran  language. 
This  enabled  the  program  to  be  written  in  essentially  the  same 
way  as  the  algebraic  solution,  if  the  letters  were  written  in 
stepwise  fashion. 


Nomenclature  for  I,  B.M.  7074  Fortran  Program 


X  =  X 

for  fixed  and  rotational  coordinates 

Y  =  Y 


K  Doj  =.  s/x 

Yoo  T  =  V* 


velocity  component  on  the  X  and  Y 
axes  for  fixed  and  rotational  coordinates 


E  =  Jap/,/* 


R:  radius  of  the  flame 


ABSV  =  \l~%* TT*~  0V 

ABSR  =  \J  X*-+fA  oy  V 5* 


The  remaining  symbols  are  self-evident  from  the  program 


logic. 
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FORTRAN  PROGRAM 


C 


Compile  Run  Fortran 
Dimension  R(250) ,  £(250) 

Flame  Trajectory  Version  1 
Call  Copy  B  (2) 

1  Print  501 

501  Format  (1H  76H  Trajectory  of  Flame  and  Its 
Experimental  Determination  in  Vortex  Combustion) 

2  Read  552,  N,  OMEG,  DT,  TMAX,  MMAX 
552  Format  (110, 3F10.0,  110) 

If  (OMEG)  20  ,  20,  22 
22  Print  502,  N 

502  Format  (lHO,4HN=,  14) 

Read  551,  (R(M),  M=l,  MMAX) 

Read  551,  (E(M),  M=1 ,  MMAX) 

551  Format  (8F10.0) 

3  Print  503 

503  Format  (1HO,  3X,  1HT,  12X,  4HXDOT,  13X,  4HYDOT, 

15X,  16X,  1HY,  14X.4HABS1V,  13X.4HABSR,  14X 
3HY/X) 

4  1=0 

5  T=0. 0 
X=0. 625 
Y  =0.0 
VX=0.0 


VY  =0. 625*OMEG 
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FORTRAN  PROGRAM,  Continued 

6  T=FLOATF(l)*DT 
PT  =  1000.  0*T 

ABSV=SQRTF(VX**2+VY**2) 

ABSR=SQRTF  (X**2+Y**2) 

TANG=Y/X 

7  Print  505,  PT,  VX,  VY,  X,  Y,  ABSV,  ABSR,  TANG 

505  Format  (1HO,  LX,  F4.  2,  3X,  7E17.  8) 

8  J=1 

9  SDVX=0. 0 
SDVY  =0.  0 
SDX=0, 0 
SDY  =0.  0 
B  =  1.0 
C=0.0 

10  TT  =T+C*DT 

VJC+VX+C*DVJX 
VJY  +VY  +C*DV  JY 
XJX=X+C*DXJX 
XJY=Y+C*DXJY 
21  M-TT 

FM  =F  LOAT  F(M) 

Ml  =M+1 
M2=M+2 

RIJ  =R(M1)  +  ((R(M2)  -R(M1))*(TT=FM))*  1000.0 
El J  =E (M 1)  +  (( E(M2)  -E(M1))*TT  -FM))*  1000. 0 
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FORTRAN  PROGRAM,  Continued 

RE=EIJ*SQRTF((VJX+XJY*OMEG)**2+(VJY-XJX* 

OMEG)**2) 

200  IF(RE-0.  1)100,  101,201 

201  IF(RE-l.O)  101,  102,202 

202  IF(RE- 10. 0)102,  103,  203 

203  IF(RE-1000. 0)103,  104,  204 

204  IF(RE-20000  0.  0)104,  105,20  5 

206  IF(RE- 1000000.  0)106,  107,  107 

100  DC  =  100. 0 
Go  to  1 1 

101  CD-10. 0*RE**(-0.  782) 

Go  to  1 1 

102  CD  =  10.0*RE**(-0.  552) 

Go  to  1 1 

103  CD=0.  557#RE#*^0.  239) 

Go  to  1 1 

104  CD=RE#*0.  0634/0. 647 
Go  to  1 1 

105  CD=25.4*RE**(-0.  238) 

Go  to  1 1 

106  CD=0.  02*RE**0.  212 
Go  to  1 1 

107  CD=0. 0 

11  DVJX=(-3.0*CD/(3. 14159*RIJ)*SQRTF((VJX+ 

XJY  *OMEG)  *#  2+ ( V  JY  -XJX*OMEG)**2)*(VJX+ 
XJY*OMEG**2)*DT 


no 


FORTRAN  PROGRAM,  Continued 

DVJY  =  (-3.0*CD/(3. 14159*RIJ)*SQRTF((VJX+ 
XJY*OMEG)**2+(VJY-XJX*OMEG)**2)*(VJX- 
XJX*OMEG)  -3.  0*XJY*OMEG**2)*DT 

DXJX=VJX*DT+0.  5*DVJX*DT**2 

DXJY  =VJY*DT+0.  5#DVJY#DT**2 

12  SDVX=SDVX+B*DVJX 
SDVY  =SDVY  +B*DVJY 
SDX=SDX+B*DXJX 
SDY  =SDY+B*DXJY 

13  Go  to  (14,  16,  15,  17),  J 

14  B=2. 0 
C=0.  5 
Go  To  16 

15  B  =  1 . 0 
C  =  1.0 

16  J=J+1 
Go  To  10 

17  VX=VX+0. 16666667*SDVX 
VY=VY+0. 16666667*SDVY 
X=X+0. 16666667*SDX 

18  TEST=TT+DT-TMAX 
IF(TEST)  19,  2,  2 

19  1=1+1 
Go  to  6 

20  Call  Copy  E 

Stop 

End 


Ill 


The  basic  program  given  here  was  used  for  the  solution 
of  Equations  2.9,  2.  10  and  2.  16  for  the  flame  trajectories 
recorded  in  Fig.  4.  16,  4.  17,  5.  4.  The  only  change  necessary 
in  the  program  was  in  statement  11,  which  required  a  change 
in  mathematical  operations  depending  upon  which  equation  was 
being  solved. 

Following  is  a  tabulation  of  the  constants  used  on  the  10 


data  cards  required  in  the  previous  Fortran  Program.  These 


data  cards  were  the  same  for  each  of  the  three  equations  solved. 

90094. 20000000.000100000. 02200000  23 

0.000000000. 300000000.440000000. 530000000 
. 6000000000. 655000000. 7080000000. 76000000 

R  0.808000000.852000000. 894000000. 930000000 

. 970000001. 000000001. 034000001.06000000 

1.090000001. 120000001. 150000001. 17600000 
1. 200000001. 220000001. 24000000 

0.0000000027. 500000037.000000044.  5000000 
50. 500000056. 00000061 . 00000066. 0000000 

70. 500000074.  500000077.600000080. 8000000 
83. 000000085. 000000086. 500000086. 5000000 
E 

87. 500000086. 000000083. 500000079, 0000000 
83. 000000054. 000000038. 500000 

1800188. 4000000.  0001000000. 00600000  7 

R  0.000000000.40000000.  70000000.  94000000 

1.060000001. 140000001. 20000000 

E  0.0000000034. 500000067. 500000083. 7000000 

101.40000097. 500000054.0000000 


